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20.  ABSTRACT  (Continued) 

were  0,  5,  10,  1^,  and  20  degrees;  vertical  wheel  loads  varied  from  1000  to 
ItOOO  N;  and  tire  deflections  were  0.15,  0.25,  and  0.35  of  the  undeflected  tire 
section  height.  The  clay  was  processed  and  placed  at  a soft  consistency  having 
a cone  penetration  resistance  of  290  kPa.  Relative  density  of  the  air-dry  sand 
ranged  from  loose  to  medium  with  corresponding  cone  penetration  gradients  of 
2.0  and  3.2  MPa/m.  Performance  of  the  powered  wheel  was  expressed  in  terms  of 
pull  force  coefficient,  side  force  coefficient,  required  input  torque  coeffi- 
cient, sinkage  coefficient,  and  wheel  slip  as  influenced  by  the  following  inde- 
pendent variables:  soil  strength,  tire  deflection  and  geometry,  wheel  load, 

and  wheel  turn  angle._The  findings  and  conclusions  of  this  report  are  as 
follows.  For  the  rar^  of  test  variables  considered  in  this  investigation,  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  system  of  combining  the  inde- 
pendent variables  into  dimensionless  parameters,  referred  to  as  mobility  num- 
bers, for  predicting  performance  of  powered  wheels  operating  at  a zero  wheel 
turn  angle  and  20  percent  wheel  slip  is  extended  to  include  wheel  turn  angles 
up  to  20  degrees  and  a broader  spectrum  of  v/heel  slips.  The  effect  of  increas- 
ing wheel  turn  angle Educes  the  pull  coefficient  developed  at  specific 
values  of  wheel  slip  and  mobility  number.  Over  the  wheel  slip  range  investi- 
gated, this  reduction  in  pull  coefficient  for  a specific  turn  angle  and  clay 
mobility  number  is  nearly  constant;  however,  for  sand  the  raze  of  decrease  of 
the  pull  coefficient  lessens  at  the  higher  slip  values.  For  given  values  of 
wheel  turn  angle  and  mobility  number,  the  side  force  decreased  as  increases  in 
power  (i.e.,  the  pull  made  gveater)  were  applied  to  the  wheel.  This  was  noted 
for  both  clay  and  sand  tests.  The  side  force  coefficient  for  clay  and  sand 
Increases  with  increases  in  the  wheel  turn  angle  (holding  the  mobility  number 
and  wheel  slip  constant);  increases  with  increases  in  the  mobility  number  (if 
the  wheel  turn  angle  and  wheel  slip  is  held  constant);  and  decreases  with  in- 
creases in  wheel  slip  (for  constant  values  of  mobility  number  and  wheel  turn 
angle).  Sinkage  coefficient  for  clay  and  sand  decreases  with  increasing  mobil- 
ity number  for  a given  turn  angle;  but  for  a given  sand  mobility  number,  the 
sinkage  coefficient  increases  with  increasing  wheel  turn  angle.  The  study  con- 
cludes by  demonstrating  the  feasibility  of  constructing  a digital  computational 
model  to  solve  nonlinear  equations  of  motion  of  a rear-wheel  drive  and  a four- 
wheel  drive  vehicle  performing  a steady-state  turn  on  soft  soil. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Sl) 

UNITS  OF  MEASUREMENT 

Units  of  measurement  used  in  this  report  can  be  converted  as  follows: 


Multiply 

By 

To  Obtain 

Metric 

(SI)  to  U.  S 

. Customary 

millimetres 

0.03937007 

inches 

centimetres 

0.3937007 

inches 

metres 

3.280839 

feet 

metres  per  second 

2.236936 

miles  (U.  S.  Statute)  per 
hour 

kilograms 

2.201*622 

pounds  (mass) 

U.  S. 

Customary  to  i 

Metric  (SI) 

inches 

25. 

millimetres 

kips  (force) 

1*1*1*8.222 

newtons 

pounds  (force) 

1*.  1*1*8222 

newtons 

pounds  (mass)  per 
foot 

cubic 

16.0181*6 

kilograms  per  cubic  metre 

CHAPTER.  I 
INTRODUCTICN 
General 

In  off-road  untracked  vehicles  the  tire  is  the  link  between  the 
ground  and  the  vehicle.  Engineering  construction  equipment,  agricultural 
and  forest  management  machinery,  and  cross-country  vehicles  having 
pneumatic  tired  tractive  devices  will  achieve  the  best  results  in  off- 
road operations  only  when  the  interaction  between  tire  and  soft  soil  is 
understood  and  quantitatively  defined  in  engineering  terms. 

A portion  of  recent  U.  S.  Army  mobility  research  was  the  develop- 
ment of  a comprehensive  analytical  model  of  ground  vehicle  systems.  In 
1971  the  existing  research  and  engineering  knowledge  of  terrain-vehicle- 
operator  interactions  was  collected,  appraised,  qnd  assimilated  into  a 
simulation  model  for  the  prediction  of  ground  vehicle  mobility  (Rula  and 
Nuttall,  1971;  U.  S.  Army  Tank-Automotive  Command,  1973).  The  second 
generation  of  this  model  was  released  Army-wide  in  1973  as  the  Army 
Mobility  Model  (AMM-75)  (Jurket,  et  al,  1975). 

Basically  the  AMM-75  consists  of  computer  modules  that  allow  for 
simulating  the  entire  vehicle  as  it  interacts  with  soil,  vegetation, 
slopes,  ditches,  obstacles,  and  other  features  of  any  geographical  area. 
The  basic  output  of  the  model  consists  of  a map  which  specifies  the 
maximum  feasible  straight-line  speed  which  the  vehicle  under  considera- 
tion might  achieve  at  any  point  in  the  terrain. 
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In  assessing  the  current  program  and  identifying  future  research 
and  development  needs  for  the  U.  S.  Amy,  the  Office  of  the  Director  of 
Defense  Research  and  Engineering  (1974)  identified  as  needed  methodology 
development  those  areas  associated  with  mobility,  agility,  and  surviv- 
ability of  combat  and  combat  support  vehicles.  Potential  need  was 
identified  for  reliable  engineering  bases  to  predict  vehicle  performance 
limits  while  maneuvering  in  off-road  terrain.  Tlie  current  computational 
modules  incorporated  in  AMC-71  and  AIM- 75  that  address  vehicle  maneuver- 
ing consist  of  simple  empirical  relations  which  do  not  address  the 
problems  in  fundamental  engineering  terms. 

One  inq)ortant  consideration  to  vehicle  maneuvering  capability  in 
off-road  terrain  derives  from  the  steering  forces  which  the  vehicle’s 
running  gear  can  generate  in  soils,  ‘fhese  forces  influence  not  only 
the  stability  of  the  vehicle,  but  also  its  power  requirements  and 
ability  to  develop  net  traction  for  slope  negotiation. 

Purpose  and  Scope 

The  principal  objectives  of  this  study  were: 

1)  Investigate  the  performance  of  single,  pneumatic  tired,  powered 
wheels  when  operating  in  the  turn  mode  one  fine-  or  coarse-grained  soils, 

2)  Formulation  of  relations  to  predict  total  side  force  developed 
by  a tire  during  a cornering  maneuver  in  soft  soil.  The  forces 
generated  result  from  the  tire  slip  angle  determined  by  the  forward 
velocity,  turn  radius,  and  wheel  steer  angle  in  combination  with  wheel 
load,  pertinent  tire  parameters,  and  characterization  of  the  soil  medium. 

3)  The  results  of  1)  and  2)  will  be  translated  into  prediction 
equations  suitable  for  incorporation  into  a digital  program  to  calculate 
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the  response  of  a four-wheeled  vehicle  executing  steady-state  flat 
turns  in  soft  soil  terrain. 

Controlled  laboratory  tests  were  conducted  with  a single  powered 
wheel  equipped  with  a 6.00-9  pneumatic  tire  and  at  turn  angles  ranging 
from  0 to  20  degrees.  Wheel  loads  and  tire  deflections  were  varied. 
Tests  were  conducted  on  a near  saturated  fine-grained  plastic  clay  of 
one  consistency  and  on  a predominantly  coarse-grained  air-dried  sand  at 
two  consistencies.  Tests  were  conducted  with  the  single-wheel  dynamom- 
eter carriage  of  the  Army  Mobility  Systems  Division,  Waterways 
Experiment  Station  (WES).  The  single  wheel  dynamometer  permits 
pneumatic-tired  wheels  to  be  tested  dynamically  at  various  controlled 
speeds  and  loads  and  under  a variety  of  consistent  and  known  soil  con- 
ditions. The  mechanical  arrangement  of  the  system  yields  measurements 
of  tire  load,  tire  deflection,  wheel  sinkage,  wheel  slip,  torque  input, 
and  net  longitudinal  and  lateral  forces. 

Current  off-road  mobility  modeling  does  not  predict  or  evaluate  a 
vehicle’s  responses  during  the  path  following  sequences  or  complex  ma- 
neuvering. The  path  following  model  should  use  as  much  of  existing 
straight-line  travel  routines  from  the  AMM  Mobility  Model  as  practical,. 
One  course  of  action  might  consist  of  an  iterative  computation  process 
of  a vehicle  moving  from  point  A to  B over  a specified  terrain.  The 
first  computation  predicts  path  performance  as  though  there  was  no 
curvature  (current  procedure),  and  with  successive  interactions  intro- 
duce curvature,  rate  of  change  in  heading  and  the  feasible  steering 
response,  and  the  required  acceleration  and  deceleration  for  obtaining 
successive  adjustments  to  speeds  along  the  actual  path.  A beginning  to- 
ward achieving  the  complete  path-following  model  is  to  develop  Initial 
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relations  between  independent  variables  of  tire,  load,  and  soil 
strength  and  dependent  performance  parameters  that  adequately  describe 
the  tire/soil  interactive  forces  generated  while  a wheeled  vehicle 
undergoes  steady-state  turning. 
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CHAPTER  II 


BACKGROUND 

General  Aspects  of  Cornering  Performance  of 
Pneumatic  Tires  in  Yielding  Soils 

General  Aspects  of  Tractive  Performance 

Tlie  vehicle  applies  forces  to  the  wheel  at  the  axle  while  the  soil 
medium  applies  forces  at  the  soil-tire  interface.  To  study  these  forces 
the  wheel  can  be  considered  as  a free  body,  disconnected  from  both  the 
vehicle  and  the  soil,  and  restored  to  equilibrium  by  forces  and  moments 
applied  to  the  axle  and  at  the  interface  area.  Figure  2.1  (after 
Schuring,  1966)  illustrates  the  possible  combination  of  torque  and 
horizontal  pull.  Vertical  force  W , horizontal  pull  P , and  torque 
M are  applied  to  the  axle,  all  of  which  are  counteracted  by  soil 
force  F . These  occurrences  can  be  illustrated  graphically  by  imagin- 
ing a wheel  with  a constant  vertical  load  W being  applied  while  it  is 
being  propelled  in  a horizontal  direction.  For  this  discussion  the 
input  torque  is  constantly  undergoing  change  beginning  with  input 
torques  that  are  opposed  to  wheel  rotation  (braked  condition)  to  positive 
values  of  input  torque  (powered  wheel) . If  these  input  torques  are 
known,  and  the  corresponding  values  of  pull  and  wheel  slip  are  measured, 
a relation  of  torque  and  pull  versus  wheel  slip  will  appear  as  qualita- 
tively shown  in  Figure  2.2.  The  various  regions  (1-5)  defined  in 
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® DRIVEN-PULLING 

© DRIVEN-FREE  ROLLING-SELF  PROPELLED  POINT 
@ DRIVEN-PULLED 

© PULLED  FREE  ROLLING  - TOWED  POINT 

® brakeo-pulled 

Figure  2.2.  Conditions  of  a Wheel  on  Yielding  Soil. 
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Figure  2.2  as  a wheel  operates  on  yielding  soils  have  been  keyed  to  the 
illustrations  on  Figure  2.1. 

The  first  quadrant  of  the  axes  shovm  in  Figure  2.2  represent  the 
operating  status  where  the  wheel  is  providing  a drawing  force  and  hence 
is  a tractive  device.  Evaluating  traction  performance  begins  by 
establishing  performance  criteria  and  determining  the  interdependency 
of  the  critical  parameters  by  making  simultaneous  measurements  in 
various  states  of  dynamic  equilibrium.  Pull,  torque,  weight  carried  by 
the  wheel,  rate  of  angular  rotation,  forward  speed,  and  wheel  slip  are 
basic  measurements  that  describe  the  state  of  dynamic  equilibrium. 

Tire  Slip  Angle 

The  lateral  forces  acting  on  pneumatic  tires  operating  on  a hard 
or  semi-rigid  roadway  have  been  frequently  investigated  and  reported. 

The  early  studies  in  the  raid-1920 's  and  early  1930 ’s  concerning  the 
mechanics  of  cornering  were  made  in  France  and  Germany.  In  the  United 
States  during  this  era  automobiles  had  solid  front  axles  and  the  center 
of  gravity  located  considerably  aft  of  the  midpoint  which  produced 
oversteer  and  instability  at  higher  speeds  (Sigel,  1966).  This 
characteristic  of  American  automobiles  plus  an  end  of  technical  exchange 
with  Germany  just  prior  to  World  War  II  produced  an  influx  of  investiga- 
tions in  this  country  during  the  1930 's  concerned  with  the  dynamics  of 
the  rolling  tire  and  relating  this  mechanical  behavior  to  the  direc- 
tional properties  of  automobiles. 

A wheel  which  is  fitted  with  a pneumatic  tire  and  which  is  con- 
strained to  a specific  plane  to  which  no  perpendicular  forces  are 
applied,  will  roll  in  a direction  coinciding  with  the  vertical  plane. 
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If,  however,  a force  is  applied  obliquely  to  the  wheel's  6ixis  then  as 
the  wheel  rolls  it  will  move  along  a path  making  an  angle  with  the  plane 
of  the  wheel.  Broulhiet  (1925)  is  credited  with  illustrating  the 
importance  of  this  angle  in  analysis  of  vehicle  handling.  Fiala  (1954) 
used  an  analytical  model  to  show  the  importance  of  the  slip  angle  con- 
cept in  generating  lateral  forces.  Segal  (1956)  combined  side  slip  with 
wheel  camber  angle  and  the  self-aligning  torque  concept  to  appraise  tire 
forces  generated  during  cornering. 

The  angle  formed  between  the  plane  of  the  wheel  and  the  instan- 
taneous velocity  vector  is  generally  designated  "slip  angle."  The  mag- 
nitude of  the  slip  angle  depends  on  many  factors.  The  most  important  of 
which  are  the  magnitude  of  the  force  applied  obliquely  to  the  plane  of 
the  wheel,  the  vertical  wheel  load,  the  inflation  pressure,  and  the 
construction  and  elastic  properties  of  the  tire  itself.  Steeds  (1960) 
states  that  the  speed  at  which  a vehicle  rolls  on  hard  surfaces  has 
little  effect  on  the  value  of  the  slip  angle. 

As  previously  implied,  an  external  force  must  be  applied  for  a 
vehicle  to  deviate  from  a straight  line.  Wlien  a tire  is  steered  across 
the  path  of  motion  a deformation  and  displacement  of  the  contact  path 
occurs  which  gives  rise  to  a side  force  (also  designated  in  the  litera- 
ture as  the  cornering  or  steering  force)  and  a moment  that  attempts  to 
realign  the  wheel  in  the  rolling  direction.  The  side  force  floes  not 
act  in  the  vertical  plane  containing  the-  axis  of  the  wheel,  but  in  a 
parallel  plane  lying  slightly  behind  the  wheel  axis,  as  illustrated  in 
Figure  2.3.  The  resulting  SF*C  acts  on  the  wheel  and  tends  to  turn  so 
its  plane  coincides  with  the  direction  of  motion;  this  is  resisted  by 
an  equal  and  opposite  couple  applied  by  steering  system  and  suspension 
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of  the  vehicle.  The  distance  C is  called  the  pneumatic  trail  and 
couple  SF*C  referred  to  as  the  trail  moment  or  self-aligning  torque. 
Self-aligning  torque  is  of  secondary  interest  in  describing  the  handling 
behavior  of  a fixed  control  vehicle,  but  is  important  which  the  loads 
within  the  steering  mechanics  are  required.  Kirch  (1973)  states  that 
the  self-aligning  torque  as  a function  of  slip  angle  is  of  little  con- 
sequence to  the  driving  behavior  of  a vehicle  in  soft  soils  because  of 
the  low  speed  encountered  with  off-road  operations. 

Driving  or  braking  a wheel  will  considerably  reduce  the  lateral 
force  obtained  at  any  given  side  slip  angle  (Ellis,  1969).  This  occurs 
because  the  additional  traction  applied  in  the  wheel  utilizes  more  of 
the  available  local  friction  which  in  turn  reduces  the  amount  available 
in  a lateral  direction. 


Related  Research  at  WES 
WES  Numeric  Prediction  System 

Tire  and  soil  modeling  is  a fundamental  part  of  the  simulation  of 
cross  country  operation  of  off-road  vehicles.  A system  was  developed 
at  WES  from  the  examination  of  a large  number  of  carefully  controlled 
laboratory  single- tire  tests  conducted  on  a saturated  plastic  clay  and 
on  two  air-dried  sands.  The  ITES  system  allows  the  prediction  of 
certain  performance  parameters  (dependent  variables)  by  combining  the 
independent  variables  'through  dimensional  reasoning  (Freitag,  1965) . 

In  the  dimensional  analysis,  the  following  independent  variables 
were  condsidered  to  influence  the  tire-soil  system  behavior 
significantly: 
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Soil:  Friction  angle 

Cohesion 
Unit  weight 
Splssltude 

Tire:  Dimensions 

Deflection 


(^) 

(c) 

ii) 

(a  measure  of  viscosity) 

(diameter,  section  width 
and  hel^t) 

(5)  as  measured  on  a rigid 
surface 


System:  Load 

Translational  velocity 
Slip 

Soll-tlre  friction 
Gravitational  constant 

Freltag  made  use  of  the  PI  Theorem  to  establish  that  14  different  Pi 
terras  were  necessary  to  describe  the  tire-soil  system  (Bridgman,  1931; 
Duncan,  1953;  and  Langaar,  1951).  The  actual  relations  that 
exist  among  Pi  terms  cannot  be  established  analytically.  Experiments 
must  be  conducted  with  the  Pi  terms  as  controlled  variables.  The 
problem  was  somewhat  simplified  by  adopting  then  established  conven- 
tional terms  for  some  of  the  Pi  terms  and  by  inspection  (from  an  exten- 
sive background  of  judgment  and  experience),  which  lead  to  10  Pi  terms 
being  written  immediately.  The  remaining  four  Pi  terms  were  developed 
by  matrix  manipulation.  In  evaluating  the  Pi  terms  the  time  dependent 
parameters  (spissitude,  velocity)  were  disregarded  and  only  purely 
frictional  (c  = 0) , or  purely  cohesive  ((j)=  0)  soils  were  considered 
and  used  in  the  tests  that  established  performance  relations  for  use  in 
each  particular  soil  type.  With  the  type  of  the  soil  limited  to 
either  purely  cohesive  or  purely  frictional,  the  soil  is  modeled  by  cone 
penetration  values,  as  follows: 

Frictional  soils:  cone  penetration  resistance  gradient  (rate 

of  increase  of  cone  index) 

Cohesive  soils:  cone  penetration  resistance 
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The  independent  variables  were  combined  by  dimensional  analysis 
into  so-called  mobility  numbers  of  the  two  soil  types  tested: 

1)  Clay  mobility  number  : 


- 1 (£\^  ^ 1 
c " W \h/  T . b 

2d 


and 

2)  Sand  mobility  number  N^,  : 


_ G(bd)^^^  . 6 
s ~ W h 


(2.1) 


(2.2) 


where 

C = Average  cone  penetration  resistance  of  the  0-  to  15-centimetre 
soil  layer  as  measured  with  WES  standard  cone  penetrometer. 

G = Average  cone  penetration  resistance  gradient  of  the  0-  to 
15-centimetre  soil  layer  as  measured  with  the  IffiS  standard 
cone  penetrometer. 

b = Unloaded  tire  section  width. 

d = Unloaded  tire  diameter. 

W Vertical  load  applied  to  the  tire  through  the  axle. 

6 = Difference  between  unloaded  and  loaded  tire  section  heights, 
h = Unloaded  tire  section  height. 

Relations  were  established  between  each  of  the  following  perform- 
ance parameters  and  the  sand  and  clay  mobility  numbers,  respectively: 
pull  coefficient  P/W  , torque  coefficient  M/Wl  , and  sinkage 


^The  dimensional  term  — did  not  appear  in  Freitag's  work  but 

1 + — 

2d 

was  added  later  by  Tumage  (1972)  to  enable  the  total  collapse  of 
additional  laboratory  data  obtained  with  tires  having  large  b/d  ratios 
(e.g.  ‘terra  tires  whose  width,  b -may  be  equal'to  the  diameter  d ). 
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coefficient  z/d  , all  at  20  percent  wheel  slip,  and  towed  force 

coefficient  P^/W 

where 

P ■ Net  pull. 

M = Torque  input  to  the  axle. 

r = Effecto-ve  tire  radius  in  the  soil. 

a 

z = Sinkage. 

These  relations  (Figure  2.4  for  clay  and  Figure  2.5  for  sand)  describe 
the  performance  at  the  towed  point  and  at  20  percent  slip.  The  pull 
generated  at  20  percent  slip  is  generally  referred  to  as  the  "maximum 
pull."  Actually  higher  pull  values  may  occur  at  other  slip  values  but 
greater  amounts  of  input  power  are  required  and  the  trade-off  between 
pull  and  required  torque  has  been  shown  to  be  optimized  at  a wheel  slip 
of  20  percent  (Freitag,  1965). 

More  recent  developments  include  expanding  the  dimensionless  single 
tire  relations  to  permit  prediction  of  tire  performance  over  a broader 
slip  range.  Smith  (1976)  performed  a thorough  reanalysis  of  the  basic 
laboratory  data  obtained  during  the  1960 *s  and  developed  relations 
describing  critical  performance  parameters  over  the  pull  slip  range  for 
single  tires  operating  in  clay.  Smith  began  by  altering  the  clay 
numeric  so  that  rigid  wheels  could  be  considered  (l.e.  “ 0 for  rigid 

wheels  and  Equation  2.1  and  2.2  would  be  equal  to  zero)  within  the 
same  framework  as  pneumatic  tires.  Smith  suggested  that  the  clay 
mobility  number  (herein  designated  as  Nj_  ) take  the  form: 
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Figure  2.5.  Hulution  of  Performance  Coefficients  to  Sand  Mobility  Humber  (after  Turnage,  1972.  Figure  3) 


Smith  found  for  powered  wheels  with  constant  test  conditions  (i.e. 
a constant  value)  a linear  relation  existed  between  the  pull 


M 

coefficient  and  the  input  torque  coefficient  (Figure  2.6a).  is 

a 

p 

defined  as  the  input  torque  at  the  self-propelled  point  (i.e.  — = 0) 

corresponding  to  the  slip  at  self-propelled,  S as  shown  in 

P 

Figure  2.6b.  By  Inspection  of  Figure  2.6a  the  relation  between  ^ and 


M 

Wr 


is  seen  to  be 


T>  rf  - M „ 

I . K.J, 

W Wr 


(2.4) 


where 


For  most  conventional  pneumatic  tires  K.J,  varies  from  about  .88  to  .97, 
A rslation  was  found  to  exist  between  the  slip  at  the  self-propelled 
point,  Sgp  and  the  clay  numeric 


0.005 


(2.5) 


As  becomes  very  large  approaches  0.005.  For  example  if  C 

became  infinitely  large  (rigid  pavement)  the  self-propelled  point  would 
be  approximately  0.5  percent  slip.  If  is  related  to  , it 

M 

would  not  seem  unreasonable  from  noting  Figure  2.6b  that  is  like- 

wise  related  to  N*  as 
c 


M 

_§£  = 
Wr 

a 


12 


(»;) 


+ 0.007 


(2.6) 
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Figure  2.6a  also  suggestes  that  the  loci  of  equal  slips  exists  as 
various  test  conditions  (i.e.  changes  in  ) occur.  Smith  found  the 
following  relation  to  best  fit  the  existing  laboratory  data. 


1°  / r- 

Equation  2.4  through  2.7  permit  ^ to  be  computed  from  given  test 
conditions  (and  likewise  the  required  torque  input)  for  a given 

slip. 

'fhe  relation  between  the  revised  clay  numeric  and  the  towed 
coefficient  — was  found  to  be 


2=0.007  (2.8) 

M 

which  is  incidentally  numerically  equal  to  , however,  the  physical 

^ a 

relation  between  the  two  dependent  variables  is  not  apparent. 

Melzer  (1973)  developed  a prediction  system  for  the  pull  and  the 
power  required  of  a single  wheel  operating  in  sand  of  various  densities. 
The  primary  results  of  earlier  studies  was  a system  for  predicting 
maximum  pull  (system  output),  torque  (system  input)  necessary  to 
develop  the  maximum  pull,  and  towed  force  (zero  torque)  for  pneumatic 
tires  operating  on  soil.  Melzer *s  relations  are  limited  because: 

1)  only  a representative  portion  of  the  sand  data  was  selected  for 
inclusion  in  the  development;  and  2)  the  slip  range  was  limited  from 
the  slip  at  the  towed  point  to  that  of  20  percent.  Only  a portion  of 
the  available  sand  data  was  reanalyzed  because  the  new  data  had  to  be 
compiled  for  values  of  pull  and  torque  at  intermediate  values  of  slip 
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between  die  self-propelled  point  and  20  percent  slip  - a time  consuming 
and  expensive  task  because  most  of  the  information  had  to  be  "hand-read" 
from  oscillographs. 

Melzer's  sunmary  relations  consisted  of  a nomograph  plot  of  power 
2 

coefficient  as  a function  of  the  pull  coefficient,  slip,  sand  mobility 
number,  and  slope  angle  of  the  soil  surface  that  the  wheel  must 
negotiate  (Figure  2.7).  Using  the  basic  data  used  to  develop 
Figure  2.7  more  fundamental  plots  of  the  output  pull  coefficient 
(Figure  2.8a)  and  input  torque  coefficient  (Figure  2.8b)  were 
constructed. 

Using  modified  Honeywell  (1971)  computer  library  routines  for  curve 
fitting  the  following  relations  were  formulated  from  the  curves  illus- 
trated in  Figure  28.  For  the  pull  coefficient: 


P ^ ^P 

“ ''s  - "V  ®P 


(2.9) 


where 


0.69 


10.8 

2.23 

„l/3 


0.01 

s 


and  the  input  torque  coefficient 


M 


where 

Ajj  * 0.66 


(2.10) 


Power  coefficient  PN  * 


M 

Wr 


a 


1 - s 
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WES  Tovmed  Turn  Tire  Test  Programs 

Green.  (1971)  reported  the  results  of  a series  of  towed  tests  con- 
ducted with  a Lunar  Roving  Vehicle  (LRV)  wheel  operating  on  a lunar  soil 
simulant  of  crushed  basalt.  Tlie  turn  angle  of  the  wheel  was  varied 
obliquely  to  the  direction  of  travel  from  -5  to  +90  degrees;  in  tlie 
latter  case,  the  plane  of  the  wheel  was  perpendicular  to  the  travel 
direction.  The  applied  wheel  loads  were  extremely  light  and  the  for- 
ward velocity  was  varied  from  1.07  to  3.05  metres  per  second. 

Side-force  coefficient  S/W  increased  with  increasing  turn  angle 
o to  a value  of  about  1.2  at  a = 90  degrees.  Speed  had  an  effect  on 
side  force  with  slightly  higher  values  of  S/W  generally  corresponding 
to  the  higher  wheel  speed. 

Unpowered  single  wheel  tests  were  conducted  at  the  WES  in  1973  to 
define  relations  between  side  force  and  turn  angle  for  two  pneumatic 
tires  common  to  the  landing  gear  structure  of  military  transport  aircraft 
capable  of  operating  on  unprepared  soil  runways  (Krick,  1975).  Tests 
were  conducted  on  mortar  sand  and  Vicksburg  buckshot  clay  of  various 
soil  strength.  Melzer  (1976)  agumented  those  tests  by  incorporating  a 
third  tire  of  differing  dimensions,  expanded  test  conditions,  and  a 
third  soil,  Yuma  sand.  All  together  Melzer  reported  on  99  one-pass 
unpowered  single  wheel  tests  conducted  in  the  laboratory  on  the  one  clay 
and  the  two  sands  with  8,50-10,  7.00-6,  and  6.00-9  tires;  turn  angles 
were  0,  5,  10,  15,  and  20  degrees.  Wheel  loads  were  varied  from  about 
1000  to  7000  N;  tire  deflections  were  0.15,  0.25,  0.35,  0.40  of  the 


23 


undeflected  tire  section  heights.  Clay  cone  penetration  resistance 
ranged  between  255  and  540  kPa.  The  air-dry  sands  had  cone  penetration 
resistance  gradients  ranging  from  0.7  to  4.6  MPa/ra. 

In  generalities  Melzer  found  that  the  side-force  coefficient  S/W 
for  clay  and  sand  increased  as  the  mobility  number  becaitu  larger  for  a 
given  turn  angle  a ; conversely,  S/W  increased  with  a when  the 
mobility  number  was  held  constant.  Also  for  both  clay  and  sand  tests 
the  towed  force  coefficient  P^/W  did  not  show  a well  defined  depen- 
dency on  the  turn  angle  a . With  consideration  of  these  observations 
Melzer 's  tabulated  data  was  used  to  formulate  relations  (between  forces 
acting  on  the  wheel  and  turn  angle  for  varied  soil  types  and  consisten- 
cies) consistent  with  the  wheeled  vehicle  modeling  needs  to  be  discus- 
sed in  Chapter  V.  • 

Figure  2.9  presents  the  towed  force  coefficient  clay  data  obtained 

3 

by  Melzer  as  a function  of  the  clay  mobility  number  and  the  turn 

angle  a . The  towed  force  is  that  force  acting  in  the  plane  of 

the  wheel.  The  curve  shown  on  each  plot  is  a pictorial  of  Equation  2.8, 

p 

IT  “-^+  0.007  <2. 8) 

("c) 

which  is  the  equation  for  predicting  the  towed  force  coefficient  for  a 
towed  wheel  following  a straight  path.  Figure  2.9  indicates  that  Equa- 
tion 2.8  anq)ly  predicts  the  performance  parameter  P.j,/W  independent 
of  turh  angle  a . 

Since  P^/W  is  independent  of  the  turn  angle  a , it  can  be  used 

3 

Melzer  reported  values  as  determined  from  tc'St  conditions  by 

Equation  2,lj  however  his  test  conditions  were  reported  vrfiich  peirmitted 
values  to  be  calculated  from  Equation  2.3. 
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to  normalize  S/W  so  that  the  effect  of  the  independent  variables, 
namely  those  variables  con^riaing  the  clay  mobility  number  and 

the  wheel  turn  angle  a , can  be  inspected.  Figure  2.10  illustrates  the 
relation  between  the  side  force/towed  force  coefficient  and  the  clay 
mobility  number  and  the  wheel  turn  angle.  For  the  specific  wheel  turn 
angles  and  test  condition  used  in  this  study  the  side  force  of  the  towed 
wheel  increased  nonlinearily  with  increases  of  the  wheel  turn  angle  and 
the  clay  mobility  number.  Hiese  data  cannot  be  effectively  extrapolated 
significantly  beyond  the  testing  limits.  If  all  variables  were  held 
constant  and  the  soil  strength  progressively  increased  (with  subsequent 
Increases  in  the  clay  mobility  number)  then  for  a given  turn  angle  the 
side  force  would  reach  some  maximum  value  and  then  decrease  to  a near 
constant  value  as  the  soil  approaches  a rigid  mass.  The  maximum  value 
would  occur,  for  a given  turn  angle,  when  an  optimized  condition  of  the 
two  interrelated  process  developed.  Resistance  develops  from  the  volume 
of  soil  undergoing  displacement  by  passive  action  of  the  turned  wheel 
partially  embedded  in  and  pushing  agaijist  the  soil.  Increased  soil 
volumes  are  involved  as  the  wlieel  sinks  deeper  into  the  soil  medium 
which  occurs  with  decreasing  values  of  soil  shear  strength  or  under 
larger  wheel  loads.  Conversely,  as  shear  strength  increases  sinitage 
decreases  but  the  passive  force  developed  by  the  wheel  required  to  over- 
come the  resisting  side  force  increases  per  unit  volume  of  soil  dis- 
placed because  of  the  higher  shearing  resistance  offered  by  the  soil 
mass. 

Curves  were  fitted  to  the  plotted  values  of  S/P„  versus  N'  for 

i c 

the  various  turn  angles  as  shown  on  Figure  2.10a  through  2.10d.  These 
curves  are  summarized  and  shown  as  solid  lines  in  Figure  2.10e.  The 
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following  equation  was  selected  to  analytically  represent  these 
graphical  relations; 


where 


^ m A 15.4  B 

P„  “ N'  - 7 + B 

T c 


(2.11) 


Thie  dashed  lines  shown  on  Figure  2.10c  are  the  corresponding  graphical 
representation  of  Equation  2.11  for  the  various  turn  angles  shown. 
Stability  problems  will  occur  for  Equation  2.11  for  very  small  values 
of  turn  angle  a . In  the  absence  of  towned  data  for  turn  angles  less 
than  5 degrees  it  would  be  recommended  that  Equation  2.11  be  linearly 
interpolated  by  the  ratio  of  the  turn  angle  in  question  to  a turn  angle 
of  5 degrees  and  with  coefficient  B equal  13.5,  corresponding  to  a a 
of  5 degrees. 

Prediction  of  the  forces  acting  on  an  unpowered  turn  wheel  being 
towed  in  soft  clay  can  be  made  by  using  Equation  2.8  for  determining  the 
clay  mobility  number.  Equations  2.9  and  2.10  for  determining  the  tow  and 
developed  wheel  torque,  respectively,  and  Equation  2.11  for  the  side 
force. 

Figure  2.11  presents  the  towed  force  coefficient  sand  data  reported 
by  Melzer  as  a function  of  the  sand  mobility  number  and  the  wheel 

turn  angle  a . The  continuous  curve  shown  on  each  plot  was  derived  from 
the  equation 


T „ . 0.83 

iT -0.015 


(2.12) 


This  equation  was  first  reported  by  Tumage  C1972)  and  developed  from  a 
large  data  base  of  single  wheels  equipped  with  various  pneumatic  tires 
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and  towed  through  two  sands  at  a zero  wheel  turn  angle.  Although  there 
is  somewhat  more  scatter  in  the  data  of  Figure  2.11  as  compared  to  the 
similar  presentation  of  the  clay  data,  Equation  2.12  can  be  effectively 
used  to  predict  the  force  in  the  plane  of  the  hub  for  a turned  wheel 
being  towed  through  sand. 

The  scatter  of  for  the  towed  tests  in  sand  prevented  the 

use  of  that  dependent  variable  for  normalizing  S/W  as  was  done  for  the 

clay  data.  Instead  the  relation  between  side  force  coefficient  S/W 

aiid  the  sand  mobility  number  N were  compared  and  graphically 

s 

presented  in  Figure  2.12.  Over  the  range  of  reported  by  Melzer 

(1976),  S/W  increases  with  increasing  turn  angle  a . For  a specific 
turn  angle,  S/W  increases  with  increasing  j however,  in  comparison 
with  the  corresponding  relation  obtained  for  clay  (Figure  2.10),  this 
increase  with  is  not  as  pronounced. 

A fundamental  equation  having  the  form  of  a rectangular  hyperbola 
was  selected  to  represent  the  S/W  versus  illustrated  in 

Figure  2.12.  The  derived  rectangular  hyperbola  is  moved  vertically 
within  the  plot  depending  upon  the  value  of  the  turn  angle.  Tlie  basic 
relation 

■y  = A + 0.83  - jj  55^^  (2.13) 

s 

where 

. . OTC  1.23 

A = 1,275  a 

is  shown  for  each  of  the  four  values  of  a with  broken  lines  on 
Figure  2.12.  Although  Equation  2.13  adequately  describes  the  experi- 
mental data  for  relating  S/W  versus  , values  of  S/W  will  result 
for  a = 0 thus  predicting  a physical  anomaly.  In  the  absence  of  towed 
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data  for  turn  angles  leas  than  5 degrees  it  would  be  recommended  that 
Equation  2.13  be  linearly  Interpolated  by  the  ratio  of  the  turn  angle  in 
question  to  a turn  angle  of  5 degrees,  thus  coefficient  A equals  0.06. 

Other  Pertinent  Studies 

Scliwangliart  (1968)  reported  on  the  results  of  tests  conducted  with 
towed  single  wheels  equipped  with  agricultural  machinery  tires  and  front 
tires  of  tractors  of  various  sizes.  Measurements  of  towing  force,  lat- 
eral force,  slnkage,  and  wheel  slip  were  determined  for  each  test  condi- 
tion. The  single  soil  used  for  this  study  was  reported  as  a sandy  clay 
processed  in  a fixed  test  bin  to  a moisture  content  of  14.5  percent  with 
the  angle  of  internal  friction  ^ varying  from  30  to  36  degrees  and 
cohesion  of  near  zero.  Figure  2.13  exemplifies  Schwanghart*s  finding  by 
illustrating  the  results  from  one  test  tire  Inflated  to  1 atmosphere 

pressure  and  tested  at  various  turn  angles  up  to  28  degrees  and  wheel 

L 

loads  between  1000  and  4000  N . Schwan^art  noted  that  the  towing  force 
in  the  plane  of  the  wheel  is  nearly  independent  of  the  wheel  turn  angle 
a (Figure  2.13a)  when  tested  at  a particular  wheel  load  and  up  to  a 
turn  angle  of  about  20  degrees.  The  side  force  (Figure  2.13b)  in  the 
plane  of  the  wheel  and  wheel  sinkage  (Figure  2.13c)  increased  with  wheel 
turn  angle  at  a certain  load.  Schwanghart  results  are  completely  com- 
patible with  the  program  on  towed  tires  in  the  turned  mode  reported  by 
llelzer  (1976)  and  discussed  earlier. 

Taylor  and  Birtwistle  (1966)  investigated  the  most  effective  wheel 


4 

Figure  2.13  is  reproduced  almost  directly  from  Schwanghart  (1968) 
and  force  units  are  designated  as  kp  where  1 kp  w 1 kg  and  1 kg  is 
equal  to  9.81  N. 
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angular  setting  of  the  disc  plough  supporting  wheels  which  would  provide 
for  a maximum  ratio  of  lateral  force  to  drawbar  force  (i.e.,  side  force 
to  counteract  side  forces  generated  by  the  plough  discs  with  minimum 
amount  of  forward  motion  resistance).  The  reported  coordinate  system 
placed  the  drawbar  force  parallel  to  the  direction  of  forward  motion  and 
side  force  perpendicular  to  the  direction  of  forward  motion.  Towed 
tests  were  performed  with  a wheel  equipped  with  7.50-16,  6-ply  tire  of 
three  tread  configuration  and  with  inflation  pressures  either  40  or 
70  psi.  Two  loads  of  1000  or  1500-pounds  were  selected  for  variation  of 
the  vertical  wheel  load.  Two  soils  were  selected  for  testing.  One  was 
a sandy  clay  of  very  .’’.ow  plasticity  (LL  =21,  PI  = 8)  and  having  an 
averaged  in-place  moisture  content  of  19  percent.  The  other  soil  was  a 
moderately  plastic  clay  (LL  = 46,  PI  = 20)  and  with  an  average  placement 
moisture  content  of  33  percent.  For  most  test  conditions  the  test 
soils  were  processed  to  two  general  conditions,  a loose  and  slightly 
compacted  consistency  so  as  to  simulate  the  wheels  of  a multi-disc 
plough  which  generally  have  two  wheels  run  the  furrows  (loose  soil  condi- 
tion) and  the  other(s)  on  the  unploughed  ground  (compacted  soil). 

Taylor  and  Birtwistle  found  that  for  the  agricultural  tires  towed 
through  the  two  soils  that  wheel  camber  and  tire  tread  pattern  had  very 
little  effect  on  the  drawbar  pull.  In  terms  of  measured  side  forces 
perpendicular  to  the  direction  of  travel,  Taylor  and  Birtwistle  results 
showed  that  the  side  force  increased  with  increases  in  wheel  turn  angle 
(all  other  variables  being  held  constant)  up  to  an  angle  of  about 
12  degrees  and  then  decreased  slightly  with  further  increases  of  the 
wheel  turn  angle.  The  effect  of  wheel  camber  (0  to  15  degrees)  was 
somewhat  linear  for  a given  wheel  loading  and  at  a specific  wheel  turn 


34 


angle  the  side  force  was  Increased  by  8 to  10  percent  for  each  five 
degree  increment  of  wheel  camber.  Magnitude  of  side  forces  was  found 
to  separate  according  to  tire  tread  pattern  for  tests  conducted  on  the 
sandy  clay  soil;  however,  this  separation  was  not  apparent  for  tests 
performed  with  the  lean  clay  test  soil. 

Krick  (1973)  reported  the  results  of  tests  conducted  with  single 
wheels  equipped  with  agricultural  tractor  and  cross  coujitry  vehicle 
tires  mounted  within  a six-degree-of-freedom  dynamometer  system.  The 
results  reported  were  conducted  within  a soil  bin  containing  a sandy 
loam  soil  moderately  compacted  to  resemble  tractor  traffic  during 
cultivation  and  at  a moisture  content  of  14.5  percent;  reference  to  a 
measure  of  soil  strength  or  its  physical  properties  was  not  provided. 

Performance  was  expressed  in  terms  of  side  forces  acting  per- 
pendicular to  the  plane  of  the  wheel  and  tractive  forces  in  the  plane 
of  the  wheel.  From  the  experimental  data  Krick  developed  relations 
between  tractive  and  side  forces  as  a function  of  side  slip  angle 
and  wheel  slip.  Example  plots  shown  in  Figure  2.14  are  typical  of  Krick 
results  developed  with  a 7.50-18  tire  at  a pressure  of  1 atmosphere 
with  a wheel  load  of  530  kp.  These  results  indicated  that  for  the 
given  test  conditions  (to  include  wheel  turn  angle)  the  side  force 
decreases  eind  the  tractive  force  increases  with  increasing  wheel  slip 
(Figure  2.14a  and  2.14b).  Krick  depicted  the  relations  between  trac- 
tive force,  side  force,  wheel  turn  angle  (designated  as  side  slip  angle 
by  Krick),  and  wheel  slip  with  "characteristics  graphs"  of  the  type  shown 
in  Figure  2.14c.  In  order  to  insure  linear  equation  of  steady  state 
motion,  Krick  linearized  the  characteristics  curves  (Figure  2.15b) 
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RACTIVE  FORCE,  KP 
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a.  LINEARIZED  TRACTIVE  FORCE  VERSUS 
SLIP  ANGLE 


b.  LINEARIZED  SIDE  FORCE  VERSUS 
TRACTIVE  FORCE  RELATIONS 


Figure  2.15»  Linearized  Turn  Tire  Test  Results. 

Wheel  Load  = 530  kips.  Tire  - 
7. 50-18  AS  Tire  with  1 atm  (98  kPa). 
Soil  = Sandy  Loam  at  ik  Percent 
Moisture  Content.  (After  Krick, 

1973). 
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and  assumed  a linear  dependence  tractive  force  and  side  slip  angle 


for  given  value  of  wheel  slips  (Figure  2.15a). 
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LABORATORY  EQUIPMENT,  TEST  PROCEDURES , 

AND  TESTING  PROGRAM 

Introduction 

This  chapter  describes  the  soils,  equipment,  and  procedures  used 
to  carry  out  the  research  on  powered  pneumatic  tires  operating  in  soft 
soils  in  a turned  mode.  Pertinent  engineering  properties  of  the  two 
soils  are  presented  and  discussed.  The  section  on  equipment  describes 
the  pneumatic  tire,  dynamometer,  cone  penetrometer,  and  related  instru- 
mentation used  to  measure  the  forces  and  important  physical  parameters 
while  conducting  turned  tire  tests.  Preparation  of  soil  cars  and 
testing  procedures  is  also  described. 

Materials 

The  entire  research  program  was  carried  out  on  two  soils  that 
represent  the  limits  of  the  soil-type  spectrum;  a near  saturated 
purely  cohesive  soil  (<|)  = 0)  an.d  an  air-dried  cohesionless  sand  (c  = 0) . 

Clajr 

One  of  the  two  soils  tested  was  a cohesive,  alluvial  clay  obtained 
from  floodplain  deposits  of  the  Mississippi  River  near  Vicksburg, 
Mississippi  and  is  locally  referred  to  as  Vicksburg  buckshot.  This 
material  is  classified  as  plastic  clay  (CH)  according  to  the  Unified 
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Soil  ClasBlflcatloa  System.  The  coaaisteiicy  data,  together  with  the 
particle  size  distribution  curve  are  recorded  in  Figure  3.1.  Compaction 
characteristics  for  Vicksburg  buckshot  are  shown  in  Figure  3.2.  The 
compaction  test  was  performed  using  standard  compaction  effort  according 
to  procedures  given  in  ASTM  Standard  D-698-70,  method  A (ASTM,  1975). 

As  indicated,  the  particular  buckshot  tested  has  an  optimum  moisture 
content  of  21.4  percent  corresponding  in  a maximum  density  of  99.1  pounds 
per  cubic  foot. 

Sand 


The  sand  used  in  the  laboratory  tests  was  taken  from  an  active  dune 
area  near  Yuma,  Arizona.  Figure  3.3  shows  the  gradation  and  index 
properties  of  this  soil,  which  is  imiformly  graded,  subangular,  and 
classified  as  SP-SM  in  accordance  with  the  Unified  Soil  Classification 
System. 

Preparation 

The  soils  were  prepared  in  movable  soil  bins  (Figure  3.4)  that  are 
0.8  metres  deep,  1.6  metres  wide,  and  long  enough  to  accommodate  test 
lanes  16  metres  long.  The  procedures  used  tc  prepare  clay  and  sand  test 
bins  with  the  desired  consistencies  and  relative  density,  respectively, 
are  briefly  described  in  subsequent  paragraphs  with  detailed  narrative 
having  been  previously  made  by  McRae  et.al.  (1965). 

Soil  Bin  Preparation,  ^oil  preparation  began  by  drying  the  soil 
to  a uniform  low  water  content  of  about  nine  percent.  Lumps  within  the 
dried  soil  were  then  reduced  by  mechanical  crushing  to  a maximum  1/4- 
Inch  size.  The  soil  was  next  blended  with  the  desired  amount  of  water 
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INDEX  PROPERTIES  OF  YUMA  SOIL 


PROPERTY 


VALUE 


SPECIFIC  GRAVITY  OF  SOLIDS 

COEFFICIENT  OF  UNIFORMITY 
CU  - Deo/D)o 

CURVATURE  COEFFICIENT 
CC  = (D3o)^  (Djo  ^ Dgo) 

D)o 

MAXIMUM  DENSITY 
DRY  UNIT  WEIGHT 
VOID  RATIO 

MINIMUM  DENSITY 
DRY  UNIT  WEIGHT 
VOID  RATIO 


i:,65 

1,50 

1,26 


0.08  MM 


104,3  LB/  CU  FT 
0.587 


87.4  LB/  CU  FT 
0,893 


Figure  3*3»  Grain  Size  Distribution  and 
Index  Values  of  Yuma  Sand. 


in  a pugmill  (Figure  3.5)  of  the  type  used  in  brickmaking  plants.  Dif- 
ferent degrees  of  soil  strength  can  be  achieved  by  preparing  the  soil  at 
different  moisture  contents.  Tlie  prepared  soil  issues  continuously  from 
the  end  of  the  pugmill  into  the  soil  bins,  which  are  propelled  slowly 
back  and  forth  by  a forklift  truck.  When  sufficient  soil  has  been 
deposited  into  the  bin  to  produce  an  approximately  six  inch  compacted 
layer,  a lumberyard  straddle  truck  equipped  with  a heavy  pneumatic 
tired  roller  (Figure  3.6)  provides  the  compacting  effort.  For  this  test 
program  sufficient  compaction  was  applied  to  achieve  a desired  compacted 
dry  density  of  88.6  pounds  per  cubic  foot  corresponding  to  a 33.0  per- 
cent moisture  content  and  a degree  of  saturation  of  99  percent.  Addi- 
tional layers  were  then  added  and  compacted  until  the  soil  bin  was 
filled.  Finally  the  surface  was  leveled  by  using  a grader  blade 
attached  to  the  straddle  truck.  After  construction  was  completed, 
measurements  were  made  to  determine  whether  the  desired  soil  conditions 
had  been  achieved. 

Previous  experience  has  fo’.md  that  the  clay  test  bins  can  be  re- 
constituted several  times  after  being  subjected  to  tire  tests  by  filling 
the  ruts  left  by  the  test  tire  and  recorapacting  with  a pneumatic  tired 
iller.  Before  each  test  five  cone  penetrations  were  made  along  the 
traverse  of  the  test  cars  to  ensure  an  acceptable  and  uniform  con- 
sistency had  been  acliieved. 

Sand  Car  Preparation.  Uniform  deposits  of  air-dried  sand  is 
achieved  by  allowing  the  soil  to  fall  through  a 1/4-lnch  mesh  screen  in 
uniform  layers  until  the  test  bins  were  filled.  The  procedure  is  il- 
lustrated in  Figure  3.7.  The  desired  test  density  is  achieved  with  the 
use  of  a small  vibrating  skid  unit  able  to  deliver  a dynamic  force  of 
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Figure  3.7*  Front-end  Loader  Spreading  Sand  on  Screen  During  Filling  of  Soil  Bins. 


1800  pounds  at  a rate  of  3600  blows  per  minute.  The  specific  density 
was  achieved  by  controlling  the  speed  of  travel  and  the  number  of  stops 
of  the  vibrator  over  the  sand  surface. 

The  objective  of  soil  processing  of  a test  bin  is  to  prepare  uni- 
form test  sections  in  which  the  Increase  in  strength  (as  determined 
with  a cone  penetrometer)  with  depth  is  approximately  linear  to  a depth 
at  least  as  great  as  the  width  of  the  test  tire.  Generally  five  cone 
penetration  determinations  were  made  along  the  test  cars  traverse  prior 
to  a tire  test.  After  a tire  test  was  performed  the  test  section  was 
rehabilitated  by  scarifying  (Figure  3.8)  to  a depth  of  150  millimetres 
and  revlbrating  until  the  desired  consistency  was  achieved. 

Test  Equipment 


Tire 


A 6.00,  4-PR  trailer  tire  buffed  free  of  tread  was  used  during 
the  testing  program.  Pertinent  tire  data  are  listed  in  Table  I.  The 
selection  of  this  tire  was  somewhat  dictated  by  the  dimensions  of  the 
modified  carriage  system  used  with  the  existing  dynamometer  system 
that  will  be  described  in  the  following  paragraphs. 

Dynamometer  System 

The  dynamometer  system,  or  test  carriage,  used  in  this  study  is 

part  of  the  basic  testing  equipment  available  at  the  WES  to  investigate 

running  gears  in  single  configuration  (McRae,  1965) . The  carriage  is 

supported  by  solid  rubber-tired  rollers  on  a pair  of  overhead  rails 

aligned  over  the  soil  bins.  These  rails  are  suspended  from  cantilever 

columns  and  cross  arms  (Figure  3.9).  The  carriage  is  towed  along  the 
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Hay  Rake  Used  in  Scarifying  Sand. 


Test  Carriage  in  Position  on  Soil'  Bins-. 


rails  by  aa  electrically  drlvea  eudlesa  cable  that  Is  fastened  fore  and 
aft  to  the  carriage  passing  over  pulleys  at  the  end  of  the  track 
system.  The  speed  of  the  towing  cable,  and  thus  the  speed  of  the 
carriage,  can  be  varied  up  to  velocities  of  about  nine  metres  per 
second.  The  test  carriage  and  the  cable  can  be  shifted  transversely 
across  the  width  of  the  soli  bin. 


TABLE  I 

DATA  FOR  6.00-9,  4-PR  GOODYEAR  TIRE 
(BUFFED  SMOOTH) 


Unloaded 

Unloaded 

Unloaded 

Unloaded 

Deflec- 

Section 

Carcass 

Section 

Inflation 

tion 

Load 

Height 

Diameter 

Width 

Pressure 

% 

N 

h,  m 

d,  m 

b,  m 

kPa 

15 

1000 

0.128 

0.516 

0.159 

52.7 

15 

2000 

0.128 

0.516 

0.160 

133.5 

15 

3000 

0.128 

0.516 

0.160 

191.5 

15 

4000 

0.128 

0.516 

0.163 

281.0 

25 

1000 

0.128 

0.516 

0.159 

14.5 

25 

2000 

0.128 

0.516 

0.159 

57.2 

25 

3000 

0.128> 

0.516 

0.159 

100.0 

25 

4000 

0.128 

0.516 

0.160 

144.5 

35 

iooo 

0.128 

0.516 

0.159 

5.0 

35 

2000 

0.128 

0.516 

0.159 

28.9 

35 

3000 

0.128 

0.51b 

0.159 

56.5 

'Ihe  carriage  consists  of  a main  structure  (Figure  3.10),  which  con- 
tains the  pneumatic  load  system,  and  a lower  frame  assembly  to  which, 

I 

under  normal  circumstances  (tests  exclusively  In  strai^t  paths),  the 

test  wheel  Is  mounted  in  such  a manner  that  it  can  be  loaded  and  powered 
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loading  system 


3.10.  Overall  View  of  Dynamometer  Carriage. 


and  yet  be  free  to  move  up  and  down.  However,  for  the  research 
described  herein  the  main  carriage  system  was  modified  so  that  the 
wheels  could  be  tested  at  various  turn  angles  (Melzer,  1976).  The 
major  modification  is  an  additional  subframe  (Figure  3.11)  that  can  be 
bolted  to  the  basic  inner  frame  (Figure  3.12)  of  the  lower  frame  as- 
sembly at  the  desired  turn  angle.  Turn  angles  can  be  varied  from  0 to 
20  degrees  in  five  degree  intervals.  In  this  configuration,  the  car- 
riage can  accommodate  wheels  with  diameters  up  to  about  .65  metres  and 
with  widths  up  to  about  .22  metres^.  The  wheels  can  be  tested  either 
powered  or  towed.  In  the  latter  case,  the  chains  that  connect  the  drive 
system  with  the  wheel  axle  (Figure  3.13)  are  removed. 

The  dynamometer  system  is  instrumented  to  measure  the  following 
quantities  continuously  during  each  test:  wheel  load,  pull  of  a 

powered  wheel  or  towed  force  of  a towed  wheel  in  line  with  the  longi- 
tudinal axis  of  carriage  travel,  lateral  forces  exerted  by  the  wheel 
on  the  inner  carriage  frame  perpendicular  to  the  direction  of  travel  of 
the  carriage,  wheel  hub  movement,  carriage  velocity,  angular  velocity  of 
the  wheel,  and  applied  torque  (powered  tests). 

The  test  wheel  axle  is  rigidly  fixed  within  the  inner  frame  of  the 
lower  supporting  frame,  and  the  inner  frame  is  suspended  from  the  outer 
frame  at  the  four  comers  by  load  cells  that  are  mounted  vertically  and 
serve  as  hinges.  Tlie  hinges  allow  the  inner  frame  to  swing  longtiudi- 
nally , but  the  movement  is  opposed'  by  a load  cell  mounted  horizontally 

between -the  two  frames  to  measure  the  horizontal  force  on  the  wheel. 

This  restraint  in  wheel  diameters  is  probably  the  major  shortcom- 
ing of  the  modified  carriage;  however,  much  larger  wheels  could  not  be 
tested  because  of  the  overall  carriage  system  was  not  designed  to 
accommodate  excessive  lateral  forces. 
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3.13.  Subframe  Mounted  in  Place  and  Wheel  Drive  System. 


In  addition,  two  load  cells  are  Installed  parallel  to  the  front  and 
rear  ends  of  the  inner  frame  and  connected  by  rods  with  the  outer  frame. 
These  load  cells  monitor  the  side  forces  exerted  by  the  wheel  on  the 
inner  frame.  Each  end  of  the  outer  frame  is  attached  to  a vertical  ball- 
spline  shaft  that  allows  the  entire  assemble  to  move  freely  in  a 
vertical  direction,  but  prevent  rotation  of  the  assembly-  in  any  plane. 

Load  is  applied  to  the  test  tire  by  means  of  pneumatic  cylinders 
mounted  between  the  upper  and  lower  frames.  Tliis  air  loading  system  is 
double-acting,  so  that  an  upward  force  can  be  used  to  permit  tests  at 
loads  less  than  the  static  weight  of  the  assemble.  The  test  carriage 
utilizes  one  pair  of  cylinders  at  the  front  and  another  pair  at  the  rear. 
The  air  storage  tanks,  which  are  visible  in  the  upper  portion  of 
Figure  3.10,  provide  a reserve  air  supply  to  compensate  for  movement  of 
the  loading  cylinders  caused  by  vertical  wheel  movement  as  it  progresses 
down  the  test  lane. 

Axial  (or  hub)  movement  is  measured  by  a potentiometer  connected 
between  the  lower  frame  assembly  (Figure  3.11)  with  the  main  carriage 
body.  Carriage  speed  was  measured  by  a tachometer.  Angular  velocity 
of  the  test  tire  was  measured  by  potentiometer  and  a tachometer  shown 
mounted  in  Figure  3.11.  Wheel  revolutions  were  monitored  by  a 
stationary  photoelectric  cell  and  a perforated  circular  disk  that 
rotates  with  the  axle,  and  carriage  position  is  indicated  by  a photo- 
electric cell  mounted  on  the  upper  frauve,  which  is  activated  by  tabs 
spaced  .10  metres  apart  on  one  of  the  overhead  rails. 

The  wheel  is  powered  by  a hydraulic  motor  driving  through  a 
specifically  constructed  mechanical  transmission  mounted  on  the  axle  of 
the  wheel.  The  wheel's  rotational  speed  can  be  regulated  at  will> 
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completely  Independent  of  the  forward  speed  of  the  test  carriage.  The 
transmission  is  restrained  from  rotating  about  the  axle  by  a connecting 
arm  in  a aeries  with  a load  cell  connected  to  the  support  frame.  During 
a test  the  input  torque  to  the  wheel  is  determined  by  recording  the 
load  ceJ.l  output  and  knowing  the  length  (moment  arm)  of  the  connecting 
members. 

Data  Recording  Equipment 

Events  measured  by  the  Inatrximents  mounted  on  the  test  carriage 
originate  as  electric  (analog)  signals  which  are  relayed  through  cables 
to  the  signal  conditioning  and  recording  equipment  (Figure  3.14).  The 
primary  recording  system  is  a FM  magnetic  tape  recorder  that  stores  the 
analog  signals  in  raw  form,  with  no  signal  conditioning,  for  further 
data  processing  (digitizing) . A secondary  recording  system  is  a 36-' 
channel,  direct-writing  oscillograph,  which  requires  signal  conditioning. 
This  latter  system  in  addition  to  providing  a backup  recording 
capability,  permits  a visual  Inspection  of  the  test  data  as  reqiiired  to 
assist  in  planning  subsequent  tests,  and  to  rapidly  appraise  teat  .. 
results.  The  accuracy  of  the  oscillograph  readings  depends  on  the 
scale  used  and  the  expertise  of  the  reader.  Only  results  obtained  from 
the  primary  recording  system  were  used  in  the  analysis  of  this  test 
program.  The  data  recorded  on  magnetic  tape  was  digitized  and  further 
processed  into  engineering  units  on  a digital  computer.  Using  appro- 
priate computer  software,  the  following  measured  parameters  were 
averaged  for  each  test;  lateral  forces  S’  , longitudinal  pull  P’  , wheel 
hub  movement,  wheel  load  W , carriage  velocity  v^  , and  translational 

a 

velocity  v of  the  wheel, 

w 
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Figure  3.1U.  Flow  Chart  for  Measurement  and  Control  Systems. 


Soil  Strength  Measurement 


Prior  to  the  conduct  of  each  tire  test,  soil  strength  was  measured 
at  five  locations  in  the  soil  cars  with  the  WES  standard  mechanical 
cone  penetrometer  (Figure  3.15).  The  WES  cone  penetrometer  was 
developed  more  than  20  years  ago  as  a device  to  obtain  an  index  of 
strength  of  surface  soils  for  trafficability  studies  and  airfield  con- 
struction. Cone  penetrometer  reading  Ci^esistance  values)  are  not  con- 
sidered basic  soil  properties  but  nontheless  a convenient  measure  of 

soil  strength,  (cone  index) . 

2 

Initially  cone  index  was  defined  as  the  average  penetration 
resistance  over  a depth  of  0 to  6 inches  in  both  cohesive  and  cohesion- 
less soils  (Green,  1964) . Later,  it  came  to  be  used  to  represent  the 
strength  of  cohesive  soils  only;  for  cohesionless  soils,  the  cone 
gradient  was  introduced,  which  is  the  rate  of  penetration  resistance 
increase  averaged  over  a depth  of  6-inches  Cfreitag,  1965).  Subsequent 
conversion  to  metric  units  results  in  the  terms  cone  penetration  and 
penetration  resistance  gradient  have  replaced  cone  index  and  cone  index 
gradient,  respectively. 

Basically,  the  instrument  consists  of  a cone  with  a base  diameter 
of  20.3  millimetres  and  an  apex  angle  of  30  degrees,  attadied  to  a shaft 
that  is  about  one  metre  long  and  has  a diameter  slightly  smaller  than 
the  cone.  A mechanized  cone  penetrometer  was  developed  for  laboratoiry 

2 

Actually  cone  index  is  a misnomer  because  the  number  is  a unit 
load  required  to  maintain  movement  of  a specifically  dimensioned  cone 
in  a soil  mass  and  actually  has  dimension  of  force  per  unit  area. 
Originally  English  units  of  pounds  per  square  inch  were  implied  but 
not  attached  to  the  number  because  the  same  size  cone  penetrometer  was 
always  used  in  related  trafficability  studies. 


Figiare  3. 15*  Mechanized  Cone  Penetrometer. 


use  In  mobility  related  research.  (McRae,  1965) . The  penetration  rate  was 
0.03  millimetres  per  second.  The  penetration  resistance  was  measured 
continuously  through  the  0-150  millimetre  depth  by  a load  cell  mounted 
at  the  top  of  the  penetrometer  shaft  and  recorded  directly  by  an  x-y 
recorder  and  simultaneously  stored  on  magnetic  tape  for  further 
processing. 

The  shear  strength  of  .soils  having  cohesive  properties  is  largely 
dependent  (disregarding  or  holding  constant  the  effect  of  previous 
stress  history,  structure  and  mineral  composition)  upon  their  density 
and  amount  of  water  present  within  the  voids.  If  a large  percentage  of 
the  void  space  of  a loose  soil  is  filled  with  air  (low  degree  of  satura- 
tion) , an  applied  load  will  result  in  compaction  of  the  soil  mass  with 
subsequent  strength  increases  (analogous  to  local  shear  failure).  If 
however,  the  voids  are  predominantly  filled  with  water  (high  degree  of 
saturation)  an  applied  load  will  be  largely  carried  by  the  pore  water 
and  a volume  change  will  occur  only  as  water  is  squeezed  from  the  mass. 
For  this  situation  the  soil  mass  would  react  to  rapidly  applied  loading 
by  yielding  when  the  cohesive  resistance  is  exceeded;  hence,  the  shearing 
strength  of  a saturated  cohesive  soil  is  independent  of  the  normal  stress 
applied.  Figure  3.16a.  Cone  penetration  of  a saturated  cohesive  clay 
will  mobilize  the  soil's  undrained  shear  strength  and  after  surface  ef- 
fects have  been  eliminated,  the  relation  between  cone  penetration  resis- 
tance and  depth  is  a unique  value,  as  illustrated  in  Figure  3.16b, 

Smith  (1966)  has  shown  that  a very  good  correlation  exists  between  cone 
penetration  resistance  as  determined  by  the  WES  penetrometer  and  cohesion 
as  determined  by  conventional  undrained-unconsolidated  triaxial  compres- 
sion tests  performed  on  clays  molded  at  sufficient  water  contents  to 

yield  degrees  of  saturation  in  excess  of  95  percent. 
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b.  PENETRATION  RESISTANCE  VS 
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c.  MOHR -COULOMB  REPRESENTATION 
OF  PURELY  COHESIONLESS  SOIL 


d.  PENETRATION  RESISTANCE  VS  DEPTH 
OF  PURELY  COHESIONLESS  SOIL 


Figure  3-l6. 


Cone  Penetration  of  Soils. 


The  available  shear  strength  of  a cohesionless  soil  is  directly  de- 
pendent on  the  applied  stress  (Figure  3.16c);  therefore,  constantly  in- 
creasing forces  must  be  applied  to  a cone  penetrometer  as  it  moves 
vertically  through  a sand  medium  (Figure  3,16d).  At  shallow  depths  this 
increased  force  is  necessary  as  the  sand's  shearing  resistance  is 
mobilized  along  the  plastic  rupture  surface  as  it  develops  while  the 
cone  moves  vertically  tnrough  a sand  medium.  Theoretically  at  some 
"critical"  depth  (the  numerical  value  depends  upon  what  method  selected, 
e.g.  Terzahi,  DeBeer,  Meyerhof t,  etc.  bearing  capacity  for  deep  founda- 
tions and  piles)  the  rupture  zone  is  fully  developed  and  penetration 
resistance  Increases  only  because  of  the  Increasing  overburden  pressure 
and  the  increase  is  therefore  much  smaller  than  the  above  the  "critical" 
depth.  Melzer  (1971)  performed  cone  penetration  tests  on  three  clean 
fine  to  medium  s^:-ids  and  found  that  the  critical  depth  using  the  WES 
cone  penetrometer  was  in  excess  of  150  millimetres  for  medium  or  dense 
sands.  This  is  below  the  depth  at  which  the  cone  penetration  resistance 
gradient  was  determined. 

The  magnitude  of  the  cone  penetration  resistance  at  any  depth  is 
determined  by  the  soil  properties.  For  soil  conditions  in  which  the 
resistance  to  penetration  is  detetmined  only  by  soil  cohesion  or  only 
frictional  properties  the  cone  penetrometer  has  been  shown  to  provide 
good  insight  to  material  properties  and  shear  strength.  However,  the 
cone  penetrometer  readings  of  soils  having  both  frictional  and  cohesion 
(most  of  the  real  world)  do  not  readily  distinguish  the  relative  effect 
of  each  consonant  of  shear  strength. 


Test  Procedures 


The  first  steps  in  the  testing  procedure  was  to  establish  zero 
positions  for  all  the  recording  tracks  and  to  record  them  both  on  the 
oscillograph  and  magnetic  tape  recorder.  The  transducer  signal  repre- 
senting each  important  variable  was  then  calibrated  to  ensure  that  the 
instrumentation  was  working  properly  and  the  calibrations  recorded. 

Before  each  test  the  soil  surface  was  leveled  and  surface  profiles 
were  taken.  Cone  penetration  resistance  was  measured  at  five  locations 
in  the  test  lane  of  the  soil  bin  before  each  test  to  check  the  uniformity 
of  the  soil  and  to  determine  whether  the  desired  soil  consistency  existed 
prior  to  testing. 

Prior  to  each  test  the  wheel  was  lowered  to  a hard-surface  platform 
adjusted  to  the  average  elevation  of  the  teat  section.  Then  the  desired 
load  was  applied  with  the  pneumatic  loading  system.  The  desired  tire 
deflection  in  percent  of  the  unloaded  section  height  (15,  25,  and  35 
percent  in  this  test  program)  was  achieved  by  measuring  the  deflected 
section  height  of  the  loaded  tire  and  adjusting  the  inflation  pressure. 

All  wheel  tests  of  this  study  were  conducted  with  a constant-slip 
technique.  The  constant-slip  rests  were  run  by  maintaining  a constant 
forward  velocity  of  the  dynamometer  system  and  a constant  angular  veloc- 
ity of  the  wheel,  by  applying  a preselected  input  torque  and  measuring 
the  pull  that  resulted.  An  unloaded  wheel  speed  of  approximately  one 
revolution  per  second  was  used  throughout  the  test  program;  the  carriage 
speed  was  adjxisted  to  obtain  the  desired  wheel  slip. 

Test  Program 


The  test  program  was  divided  into  two  parts  consisting  of  23  tests 
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performed  on  near  saturated  plastic  clay  and  49  tests  performed  on  air- 
dried  Yuma  sand. 

Constant-slip,  one-pass,  powered  tests  were  conducted  in  the 
laboratory  with  a 6.00-9,  4-PR  tire.  Performance  was  measured  in  terms 
of  pull,  side  force,  torque,  and  sinkage.  VIheel  load  was  varied  between 
1000  and  4000  N.  Tire  deflections  were  15,  25,  and  35  percent  of 
unloaded  tire  section  hei^t.  Cone  penetration  resistance,  C , was 
approximately  290  kPa  for  the  clay,  and  two  cone  penetration  resistance 
gradients  of  2.0  and  3.2  mPa/ra  in  the  Yuma  sand. 

The  tests  conducted  in  this  program  are  tabulated  in  Table  II  and 
test  results  are  summarized  in  Table  A.l,  Appendix  A,  for  clay  tests 
and  Table  A. 2,  Appendix  A,  for  sand  tests.  Figure  A.l  is  provided  as  a 
descriptor  and  key  for  the  table  headings. 


No.  of 
Tests 


4 

1 

3 

3 

3 

3 

3 

1 

1 

1 


3 

4 
6 
2 

1 

3 

1 

3 

1 

1 

1 

3 

4 

4 

4 

3 

1 

1 

3 


TABLE  II 
TEST  CONDUCTED 


Design 

Wheel 

Wlieel 

Tire 

Design 

Turn 

Load 

Deflection 

Soil 

Angle 

N 

% 

Strength 

Degrees 

Clay  Test 

kPa 

2000 

35 

290 

0 

2000 

35 

290 

5 

2000 

35 

290 

10 

2000 

35 

290 

20 

2000 

15 

290 

5 

2000 

15 

290 

15 

2000 

25 

290 

10 

1000 

25 

290 

10 

4000 

25 

290 

10 

3000 

35 

290 

5 

Sand  Tests 

mPa/m 

2000 

35 

2.0 

0 

2000 

35 

2.0 

5 

2000 

35 

2.0 

10 

2000 

35 

2.0 

20 

1000 

35 

2.0 

15 

1000 

35 

2.0 

20 

3000 

35 

2.0 

15 

2000 

15 

2.0 

10 

2000 

15 

2.0 

20 

2000 

25 

2.0 

10 

2000 

25 

2.0 

15 

2000 

35 

3.2 

5 

2000 

35 

3..2 

15 

1000 

35 

3.2 

10 

1000 

35 

3.2 

15 

1000 

35 

3.2 

20 

3000 

35 

3.2 

15 

2000 

15 

3.2 

15 

2000 

25 

3,2 

15 
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CHAPTER  IV 


ANALYSIS  OF  TEST  DATA 

Freltag  (1965)  showed  through  the  use  of  dimensional  reasoning  that 
the  important  parameters  in  predicting  the  performance  of  a powered 
wheel  operating  in  soft  soil  were  wheel  load  W ; pneumatic  tire  factors 
of  deflection  6 , diameter  d , and  width  b j and  the  strength 
characteristics  of  the  soil  as  expressed  by  cone  penetration  resistance. 
The  most  important  performance  parameters  are  input  torque  M , rim 
pull  -P  , and  sinkage  z . 

M , P , z = F (C  , W , fi  , d , b)  C4.1) 

Wheel  slip  should  also  be  considered  an  independent  variable  as  done  by 
Smith  (1975)  for  clay  and  Melzer  (1974)  for  sand.  For  a turned  powered 
wheel  an  additional  independent  variable  turn  angle  a (or  more  cor- 
rectly, effective  slip  angle)  would  be  added  and  an  additional  dependent 
variable  side  force  S would  complete  the  performance  parameters. 

Tests  in  Clay 

Performance  Parameters  of  Pull  and  Torque 

Sufficient  test  data  were  not  compiled  to  permit  the  Incorporation 
of  turn  angle  a as  an  independent  variable  into  the  clay  mobility 
uunber.  Rather,  the  effect  of  turn  angle  on  pull  and.  input  torque  was 
determined  by  developing  a con^arison  of  clay  mobility  numbers  as 
determined  from  test  conditions  and  'as  would  be  computed  from  measured 


performaace  pai^ametejrs , This  is  illustratad  iji  Figure  4.1. 

The  performance  prediction  equations  presented  in  Chapter  II  can 
be  used  to  compute  the  clay  mobility  number  if  two  of  the  three 
performance  measures  (pull,  input  torque,  and  wheel  slip)  are  known. 
Three  combinations  of  measured  test  values  are  therefore  possible  for 
computing  the  clay  mobility  number:  1)  pull  and  wheel  slip,  2)  input 

torque  and  pull,  and  3)  input  torque  and  wheel  slip.  Since  pull  and 
wheel  slip  might  be  considered  as  system  output  from  the  input  torque, 
that  combination  was  chosen  to  compute  clay  mobility  numbers  with,  the 
following  relations; 


and 


S 

sp 


+ 0.005 


C2.7) 


(2.5) 


solving  for 

N'  = 
c 

It  should  be  noted  that  using  any  two  of  the  three  measured  test 
results  of  input  torque,  pull,  and  wheel  slip  and  corresponding  rela~ 
tions,  the  computed  values  of  clay  mobility  numbers  did  not  differ 
appreciably  as  seen  in  the  tabulation  in  Appendix  B. 

Figure  4.2  illustrates  graphically  the  relation  between  the  clay 
mobility  number  computed  from  the  independent  variable  (i.e.  C , W , 
6 , d,  b)  and  the  two  chosen  dependent  variables  of  pull  and  slip  for 
each  test  condocted.  Straight  lines  have  been  fitted  po  the  plotted 
data  per  turn  angle  a . The  slope  of  these  linear  lines,  designated 


21 


.2/5 


L T 2P 

Vlog  TT 


- 0.005 


(4.2) 
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as  plotted  for  respective  turn  angles  (Figure  4.3).  A 

curve  was  fitted  to  the  data  to  permit  computation  of  intermediate 
values  of  wheel  turn  angle  a . 

N'  ^ = (1  - 2.26  • N'  (4.3) 

ADJ  c 

Hence,  given  the  independent  variables  of  tire  size  and  deflation,  soil 
strength,  and  wheel  loading  can  be  con5>uted  from  Equation  2.3  and 

with  the  turn  angle  a knoxni,  confuted  from  the  above 

relation  which  in  turn  is  used  with  Equations  2.4  through  2.7  for  com- 
puting predicted  parameters. 

Side  Forces  Developed  in  Clay 

Pull  of  a powered  wheel  operating  on  a straight  line  path  is 
proportionate  to  the  clay  mobility  number  and  desired  wheel  slip 

(determined  by  the  input  torque).  For  x^heels  in  a turned  mode  the 
wheel  turn  angle  a constitutes  an  additional  input  variable  required 
to  describe  the  pull.  It  would  appear  intuitive  that  the  horizontal 
side  force  acting  normal  to  the  hub  of  a turned  wheel  might  be  defined 
by  values  of  the  clay  mobility  number,  wheel  slip,  and  turn  angle  or 
those  same  factors  on  which  pull  is  dependent.  Hence  pull  and  side  force 
should  be  related  for  given  values  of  the  turn  angle. 

Figure  4.4  is  a plot  of  the  side  force  coefficient  S/W  versus 
the  pull  coefficient  P/W  for  a group  of  tests  having  approximately 
equal  clay  mobility  numbers  between  15  and  18.  Also  plotted  (with 

solid  symbols)  are  data  points  obtained  from  towed  wheels  in  a turned 
mode  as  reported  by  Welzer  (1976)  which  had  the  stated  r^ge  of  clay 
mobility  numbers.  Numbers  beside  the  plotted  points  indicated  percent 
wheel  slip.  Based  On  the  limited  tests  performed  linear  relations 
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ADJ 


0.087  0,175  0.262  0.349  0.436 

RADIANS 
TURN  ANGLE  « 


Figure  1|,3.  Relation  'between  and  Turn  Angle  a . 


were  assumed  between  sida-fo^ce  and  pull  coefficient  for  corresponding 
values  of  turn  angle.  These  data  would  indicate  that  for  a given  wheel 
turn  angle  and  clay  mobility  number,  the  effect  of  increasing  the  input 
power  (input  torque)  to  a wheel,  and  thereby  increasing  the  pull,  results 
in  a decrease  in  the  side  force.  The  rate  of  reducing  the  side  forces 
magnitude  increases  with  increased  values  of  the  turn  angle  a . Also 
from  the  Lowed  point  to  a slip  value  of  approximately  20  percent  this 
relationship  can  be  approximated  by  a linear  function.  With  sufficient 
data,  isobars  of  equal  wheel  slips  seemingly  could  be  constructed.  These 
data  have  that  general  form  of  the  results  reported  by  Krick  (1973)  and 
reviewed  in  Chapter  II  except  that  these  data  indicate  a more  pronounced 
increase  in  the  slope  ratio'  of  side  force  to  pull  as  the  wheel  turn 
angle  increases. 

Figure  4.5  illustrates  clay  tost  data  in  terms  of  side  force 
coefficient  versus  pull  coefficient,  where  the  wheel  turn  angle  was 
10  degrees.  Again,  towed  data  having  a turn  angle  of  10  degrees  and 
respective  values  of  clay  mobility  numbers  were  extracted  from  that 
reported  by  Melzer  (1976).  Linear  relations  have  been  placed  through 
the  data  for  wheel  slips  between  the  towed  point  and  approximately 
20  percent  slip.  These  data  indicate  that  for  a given  wheel  turn  angle 
a the  relation  of  side  force  coefficient  versus  pull  coefficient  are 
parallel  lines  displaced  higher  along  the  ordinate  as  the  clay  mobility 
number  Increases. 

It  would  seem  reasonable  to  assiime  that  as  100  percent  slip  (wheel 
in  full  spin)  is  approached  the  side  force  will  approach  zero  while  the 
pull  coefficient  will  assume  some  finite  value  approximated  by 
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SIDE  FORCE  COEFFICIENT  S/W 


Equation  2.7.  Pull  coefficients  of  0.39,  0.88,  and  1.0  correspond  to  a 
wheel  slip  of  100  percent  and  clay  mobility  numbers  of  7,  18,  and  30, 
respectively.  These  pull  coefficient  values  form  the  end  points  for 
the  three  sets  of  data  illustrated  in  Figure  4.5.  The  dashed  lines 
illustrate  a potential  path  of  S/W  versus  P/W  at  slip  values  in 
excess  of  20  percent. 

From  the  above  discussion  two  assumptions  are  made  for  justifying 
the  remaining  development  of  predicting  side  forces. 

1)  The  relation  between  S/W  and  P/W  at  any  turn  angle  is 
linear  for  wheel  slip  values  between  the  towed  point  and  20  percent 
(Figure  4.4). 

2)  For  a given  turn  angle  the  relation  between  S/W  versus 
P/W  is  linear  for  a given  clay  mobility  number  and  parallel  but 
vertically  displaced  as  the  clay  mobility  number  is  varied  (Figure  4.5). 

Slopes  of  the  three  lines  shown  in  Figure  4.4  are  plotted  against 
the  respective  turn  angle  a in  Figure  4.6.  A straight  line  was  passed 
through  the  origin  having  a slope  of  1.72. 

The  value  of  side  force  coefficients  when  the  pull  coefficients 
are  zero  were  read  from  Figures  4.4  and  4.5  and  plotted  in  Figure  4.7 
and  4.8,  respectively.  Figure  4.7  illustrates  that  the  side  force 
coefficient,  when  the  pull  is  zero  and  the  clay  mobility  number  is 
constant,  Increases  with  increasing  values  of  turn  angle  a . The 
trend  indicated  that  for  higher  values  of  turn  angle  a , the  S/W  at 
zero  pull  may  decrease.  Figure  4.8  indicates  that  if  the  turn  angle 
is  held  constant,  and  for  the  test  conditions  upon  which,  these  findings 
are  based,  S/W  at  zero  pull  increases  as  the  clay  mobility  number 
increases.  This  trend  would  not  continue  indefinitely  but  most  likely 
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SIDE  FORCE  COEFFICIENT  AT  ZERO  PULL  S/W, 


Figure 


1,7.  Side  Force  Coefficient  at  Zero  Pull  Versus 
Turn  Angle  a . Data  Points  Extrapolated 
from  Figure  k.h.  Clay  Mobility  Number 
N*  = 15  to  18  . 
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would  begin,  to  decrease  at  some  nominal  value  of  the  clay  mobility 
number  . As  approaches  an  Infinitely  large  lumber,  as  would 

occur  when  the  operating  surface  approached  a hard,  semi-rigid  medium 
(i.e.  cone  penetration  resistance  becomes  very  large),  the  side  force 
would  be  reduced  primarily  because  the  friction  properties  of  the  tire- 
surface  would  dominate  performance  rather  than  from  passive  earth  pre- 
sure  failure  and  side  friction  on  the  tire  that  must  accompany  a turned 
tire  partially  sunk  as  it  comers  in  a plastic  soil  medixm. 

The  relations  developed  from  Figures  4.7  and  4.8  permit  the  side 
force  coefficient  at  zero  pull,  per  increment  of  turn  angle,  to  be 
computed  for  any  clay  mobility  number.  For  a turn  angle  of  10  degrees: 


„ a = 10  /(n'  + 4)^ 

Wp.o 


(4.4) 


and  for  a clay  mobility  number  of  approximately  15  to  18: 


<?  = 18  , 

§ = 3.37  a - 4.24  (4.5) 

P = 0 

From  Figure  4.4  when  = 18  and  a = 10  degrees  (S/W)p  _ q is  equal 
to  0.46  which  is  used  to  normalize  when  the  two  previous  equations  are 
combined  to  give  the  side  force  coefficient  at  zero  pull  per  turn  zingle- 
and  clay  mobility  number,  or: 

S 
W 

P = 0 

The  side  force  coefficient  at  zero  pull  constitutes  the  ordinate 
intercept  with  1.72  ct  , where  a is  expressed  in  radians,  expressing 
the  slope  of  the  S/W  versus  P/W  relation.  Hence  for  a given  turn 


O.OO81) 


\ 1/2 


3.37  o - 4.24  a" 


/ 


0.46 


(4.6) 


angle  and  clay  mobility  number; 

S S 
W " W 

P = 0 


p o 

(1.72  a)  i- 


(4.7) 


Sinkage  in  Clay 

Sinkage  coefficient  z/d  (vertical  hub  movement  value  divided  by 
the  inflated  but  unloaded  wheel  diameter)  ac  a function  of  the  clay 
mobility  number  and  wheel  turn  angle  is  depicted  in  Figure  4.9.  The 
line  drawn  on  each  plot  of  Figure  4.9  represents  the  cumulative  results 
of  previous  test  programs  in  which  numerous  combinations  of  the  indepen- 
dent variables  were  considered;  however  the  wheels  were  always  at  zero 
turn  angle  and  the  wheel  slip  very  close  to  20  percent  (Tumage,  1972). 
The  data  from  this  program  arc  shown  as  plotted  points  with  the  respec- 
tive wheel  slip  noted  beside  each  netum  point.  For  a given  turn  angle, 
z/d  decreases  with  increasing  as  expected.  Insufficient  data 

prevents  meaningful  analysis  as  to  the  influence  of  turn  angle  on 
sinkage  at  a specific  value;  although  the  data  suggests  that  a is 
of  secondary  in^ortance  in  that  the  previous  described  relation  reported 
by  Turnage  (1972)  amply  depicts  results  from  this  study.  This  data  also 
suggests  that  over  the  slip  range  (0  to  20  percent)  used  in  this  test 
program,  wheel  slip  does  not  appreciably  influence  sinkage. 


For  a given  set  of  independent  variables  expressed  as  the  clay 
mobility  number  and  the  wheel  turn  angle,  forces  acting  on  a 

pneumatic  tire  are  computed  as  follows; 

The  adjusted  clay  mobility  number  accounting  for  turn  angle  o 
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(4.3) 


"abJ  -“it-  2-26 


Slip  at  the  self  •'propelled  point 


(^sp)“ 


21 


(^ADj) 


2.5 


+ 0.005 


The  input  torque  coefficient  at  the  self-propelled  point 


/m  Y 

K/  K. 


(^ADj) 


+ 0.007 


The  pull  coefficient  in  the  plane  of  the  wheel 

. 0.5  log 


(2.5) 


(2.6) 


(2.7) 


Input  torque  coefficient  is  determined  by  equating  Equations  2.6  and  2.7 

, .r . (,  .kf"' liY 

Wr 


and  finally  the  side  force  coefficient 


(f)°  ' (Dp  , „ - “)  (l)“  W-2 

Values  of  actual  laboratory  test  conditions  were  substituted  into 
the  above  equations  to  obtain  predictions  of  the  performance  coefficients 
pull,  side  force,  and  input  torque.  Table  III  provides  a tabulated 
conqiarison  of  computed  coefficients  with  those  measured  during  the 
course  of  testing.  The  comparisons  provided  in  Table  HI  were  determined 
by  taking  the  absolute  difference  of  the  corresponding  predicted  and 
measured  performance  coefficient  divided  by  the  larger  of  che  two  terms 
and  expressed  as  a percentage.  Considering  all  23  tests  conducted  in 
this  program,  74  percent  of  the  pull  coefficients,  89  percent  of  the 
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TABLE  III  (CONTINUED 


side  force  coefficients,  and  70  percent  of  the/input  torque  coefficients 
Imd  percent  differences  of  20  percent  or  less.  Figure  4.10  through  4.12 
contain  the  plotted  test  values  of  pull,  side  force,  and  input  torque 
coefficients  versus  percent  vheel  slip.  Tests  having  similar  test  con- 
ditions were  grouped  to  permit  meaningful  representation.  Superimposed 
are  lines  computed  from  prediction  Equations  2.5,  2.6,  2.7,  4.3,  and 
4.7.  Input  values  to  these  equations  were  determined  from  the  desired 
test  conditions  as  provided  in  Table  I which  differ  slightly  from  the 
individual  or  average  test  conditions  being  represented.  Effective 
work  is  not  developed  by  a powered  wheel  until  positive  pull  is  realized. 
This  program  v/as  concerned  with  the  performance  of  a po'/jered  wheel 
between  the  self-propelled  slip  and  about  20-percent  wheel  slip;  hence, 
performance  relations  are  not  shown  on  Figures  4.10  through  4.12  for 
wheel  slip  values  less  than  that  occuring  when  the  pull  is  zero. 

Sand  Tests 

Performance  Parameters  of  Pull  and  Torque 

Each  condition  for  the  powertd  turned  wheel  tests  performed  in 
Tuma  sand  (Table  A. 2,  Append?.x  A'  can  be  expressed  in  terms  of  a sand 
mobility  number  Ng  , and  wheel  slip  s , and  wheel  turn  angle  a . 
Substituting  values  of  P/W  , M/Wr  , and  s from  each  test  into 
Equation  2.9  and  2.10,  sand  mobility  numbers  were  detemined  per  test 
and  reported  in  Table  B.2,  Appendix  B.  These  computed  sand  mobility 
numbers  combine  those  independent  variables  (i.e.  test  conditions)  that 
yielded  predictions  of  performance  parameters  P/W  and  M/Wr  for  a 
powered  wheel  traveling  without  a turn  slip  angle.  With  the  complete 
absence  of  test  variation,  experimental  error,  and  difference  between 
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Influence  of  Turn  Angle  on  Performance  Parameters.  Tire  Deflection 
6/h  = 0.35  1 Design  Wheel  Load  W = 2,000  W . 


ice  of  Turn  Angle  and  Wheel  Load  on  Performance  Parameters. 
Deflection  6/h  = 0.25  > Turn  Angle  a = 10  Degrees. 


experimental  curves  and  generalized  equations,  the  Ng  computed  from 
Equations  2.9  and  2.10  would  be  identical.  Unfortunately,  such  is  not 
the  nature  of  man  or  machine.  Since  pull  and  wheel  slip  are  considered 
as  output  from  the  respective  torque  applied  to  the  wheel,  the  sand 
mobility  number  computed  from  test  values  of  P/W  and  s are  plotted 
in  Figure  4.13  against  N_  as  computed  from  actual  test  conditions  and 
grouped  accordingly  to  wheel  turn  angle.  Since  the  computed  sand 
mobility  numbers  are  derived  from  pull  coefficients  on  powered  wheels 
having  the  added  resistance  of  side  forces,  these  values  would  seemingly 
be  less  than  test  conditions  would  suggest  for  a powered  wheel  underway 
with  the  absence  of  a turn  angle.  Furthermore,  large  differences  would 
be  expected  to  occur  as  the  wheel  turn  angle  increased  and  larger  side 
forces  developed.  Hence  the  plotted  values  of  computed  versus  test  con- 
dition sand  mobility  numbers  in  Figure  4.13  should  lie  below  the  45 
degree  diagonal  lines. 

Examination  of  Figure  4.13  reveals  that  for  wheel  turn  angles  up  to 
10  degrees,  the  sand  mobility  number  determined  from  measured  pull  and 
slip  are  often  greater  than  those  Ng  expressed  by  test  conditions. 
Furthermore,  a distinct  relation  (neither  linear  like  that  noted  in 
Figure  4.2  for  the  clay  data  nor  otherwise)  for  each  wheel  turn  angle 
between  the  two  values  of  sand  mobility  numbers  is  not  apparent  although 
a general  trend  does  exist  for  larger  differences  between  the  two  sand 
mobility  numbers  being  developed  as  the  wheel  turn  angle  a increases. 

Since  the  sand  turn  tire  data  would  not  lend  itself  to  analysis  with 
the  generalized  procedure  just  discussed  a more  specific  method  was 
selected  that  involved  developing  the  pull  and  torque  data  to  investigate 
the  variation  of  the  coefficients  of  Equations  2.9  and  2.10  as  affected 
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by  the  influence  of  the  side  slip  angle.  The  procedure  began  by  placing 
smooth  curves  through  data  points  of  P/W  and  M/Wr^  versus  wheel 
slip  s.  ; tests  having  approximately  equal  sand  mobility  numbers  were 
plotted  together  for  meaningful  comparisons.  These  curves  served  as  the 
bases  for  developing  more  fundamental  relations  of  P/W  and  M/Wr^ 
versus  N_  at  certain  values  of  wheel  slip;  in  turn  these  curves  were 
measured  against  the  respective  relation  at  zero  wheel  turn  angle  to 
permit  a comparative  analysis. 

Measured  values  of  pull  and  torque  obtained  during  the  powered 
turn  tire  tests  performed  in  Yuma  sand  are  presented  in  therms  of  P/W 
and  M/Wr  versus  wheel  slip  s in  Figure  4.14  through  4.19. 
Performance  parameters  in  terms  of  P/W  and  M/Wr  from  the  various 

a 

tests  are  comparative  only  if  the  independent  variables  are  alike  or 
approximately  equal.  The  independent  variables  are  represented  in  terms 
of  the  sand  mobility  number  and  increase  from  3.5  to  26.3  in  six  unequal 
increments  as  designated  in  Figure  4.14  to  4.19,  respectively.  The 
dashed  line  drawn  through  each  figure  represent  predicted  values  of  P/W 
and  M/Wr  at  zero  wheel  turn  angle  across  the  slip  ranges  investigated 

3i 

in  this  test  program  using  Equations  2.9  and  2.10,  respectively.  The 
solid  lines  were  judged  to  best  fit  the  plotted  data  points  and  experi- 
mental trends. 

Examination  of  pull  coefficient  versus  wheel  slip  at  specified 
sand  mobility  number  and  wheel  turn  angle  indicates  that  the  pull 
increases  with  increased  values  of  slip.  The  effect  of  increasing  the 
turn  angle  is  seen  to  reduce  the  pull  at  specific  values  of  wheel  slip. 
The  generrl  slopes  of  the  P/W  versus  s plots  indicate  that  at  the 
higher  wheel  slip  values  the  rate  of  decrease  in  pull  begins  to  diminish 
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TORQUE  COEFFICIENT  M/wr^  PULL  COEFFICIENT 


a.  PULL  COEFFICIENT 


b.  TORQUE  COEFFICIENT 

Figure  4.l4.  Influence  of  Turn  Angle  on  Pull  and  Torque  Coefficients. 

Average  Cone  Penetration  Resistance  Gradient  G = 2.0  , 

Tire  Deflection  6/h  = 0.15  > Design  Wheel  Load 

VJ  = 2,000  N . Desired  Sand  Mobility  Number  = 3.5  • ' 
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TORQUE  COEFFICIENT  M/Wh^  PULL  COEFFICIENT 


LEGEND 


SYMBOL  TURN  ANGLE. a 

O 0*  TEST  DATA 


d.  PULL  COEFFICIENT 


b.  TORQUE  COEFFICIENT 


) Figure  h,l6.  Influence  of  Turn  Angle  on  Pull  and  Torque  Coefficients. 

. Average  Cone  Penetration  Besistance  Gradient  G = 3.2  , 

i Tire  Deflection  6/h  =0.25  , Design  Wheel  Load 

W = 2,000  W . Desired  Sand  Mobility  Number  =9.4  . 
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TORQUE  COEFriClENT  M/wr  PULL  COEFFICIENT 


WHEEL  SUP  S,  PERCENT 

a.  PULL  COEFFICIENT 


Figtire  4.17.  Influence  of  Turn  Angle  on  Pull  and  Torque  Coefficients 
Average  Cone  Penetration  Resistance  Gradient  G = 3.2 
Tire  Deflection  d/h  =0.35  > Design  Wheel  Load 
W = 2,000  W . Desired  Sand  Motility  Numher  =13.2  . 


TORQUE  COEFFICIENT  M/Wr^  PULL  COEFFICIENT 


TORQUE  COEFFICIENT  M/wr^  PULL  COEFFICIENT  p/ W 


Figure 


oi.  PULL  COEFFICIENT 


.19.  Influence  of  Turn  Angle  on  Pull  and  Torque  Coefficients 
Average  Cone  Penetration  Resistance  Gradient  G = 3.2 
Tire  Deflection  6/h  = 0.35  , Design  Wheel  Load 
W = 1,000  N . Desired  Sand  Mobility  =26.3  . 


for  any  turn  angle.  If  at  100  percent  slip  the  side  force  Is  reduced  to 
zero,  then  on  any  turn  angle  the  P/W  versus  s relation  will  approach 
the  respective  P/W  versus  s relations  for  a zero  turn  angle  as  the 
wheel  slip  Increases  (of  course  comparative  sand  mobility  numbers  are 
necessary) . Values  of  the  pull  coefficient  were  Interpreted  from  the 
P/W  versus  s relations  shown  In  Figures  4.14a  to  4.19a  at  wheel 
slips  of  5,  7.5,  10,  15,  and  20.  The  values  of  P/W  for  each  of  the 
particular  wheel  slips  with  the  corresponding  turn  angles  and  sand 
mobility  numbers,  were  plotted  (Figures  4.20  through  4.24)  In  the  form 
of  P/W  versus  Ng  and  curves  drawn  through  points  representing  simi- 
lar turn  angles.  Equations  of  the  form 


P 

W 


a 


A 


a 


A B 

a a 

N_  - C + B 
S a a 


(4.8) 


were  established  for  each  of  the  plotted  curves  with  the  aid  of  computer 
programming.  Notice  that  Equation  4.8  Is  the  Identical  form  to  Equa- 
tion 2.9,  which  was  established  to  predict  the  pull  coefficient  for  any 
sand  mobility  number  and  wheel  slips  up  to  20  percent.  The  final  task 
was  to  develop  relations  to  define  the  influence  of  independent  vari- 
ables of  N„  , S , and  a on  the  coefficients  A , B , and  C 
S’’  a ’ a ’ a 

of  Equation  4.8  which  differed  In  magnitude  from  the  related  coeffi- 
cients in  Equation  2.9. 

Table  IV  lists  the  values  of  coefficients  A , B , and  C from 

a a a 

Equation  4.8  for  various  wheel  slips  and  wheel  turn  angles.  Figures 

4.25  to  4.27  illustrate  how  coefficients  A , B , and  C vary  with 

a a a 

wheel  slip  and  turn  angle.  Also  shown  are  the  equations  established  to 

approximate  the  plotted  points  Into  continuous  functions.  Hence  with 

the  coefficients  defined.  Equation  4.8  to  predict  the  pull  coefficient 
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PULL  COEFFICIENT  P/W 


0.6  r- 
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SAND  MOBILITY  NUMBER'  Nj 


Figure  4.20. 


Pull  Coefficient 
V/heel  Slip  = 5 


Versus  Sand  Motility  Numter, 
Percent . 
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PULL  COEFFICIENT  P/W 


PULL  COEFFICIENT  o/W 


PULL  COEFFICIENT  P/W 


PULL  COEFFICIENT  P/W 


1 


F-gure  H.24.  Pull  Coefficient  Versus  Sand  Mobility  Number, 
Wheel  Slip  = 20  Percent. 


for  various  combinations  of 


, slip,  and  a becomes 


p 

w 


a 


A 

a 


A B 

g g 

N_  - C + B 
Sag 


where 


A = 
g 


Ap  - 1.42  g 


(0.6  + 3.1  s) 


B = 
g 


Bp  - 16  g 


0.8 


and 


C = 
g 


_ . 15  3 

S + S"“ 


(4.8) 


TABLE  IV 

VARIATION  OF  COEFFICIENT  USED 
IN  EQUATIONS  2.9  and  4.8 


Turn 


Wheel  Angle 

Slip  g Coefficients 


% 

Degree 

A 

B 

C 

5.0 

0 

0.490 

10.8 

6.05 

10 

0.142 

7.22 

6.0 

15 

0.0379 

3.36 

14.0 

20 

— 

— 

“ 

7.5 

0 

0.557 

10.8 

5.3 

10 

0.207 

5.21 

5.5 

15 

0.0907 

3.54 

10.2 

20 

0.0396 

3.46 

13.1 

10.0 

0 

0.590 

10.8 

4.80 

10 

0.285 

6.05 

5.0 

15 

0.173 

6.51 

8.3 

20 

0.0817 

3.34 

11.0 

15.0 

0 

0.623 

10.8 

4.20 

10 

0.418 

6.56 

4.5 

15 

0.290 

6.07 

5.7 

20 

0.166 

5.57 

7.7 
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TABLE  IV  (CONTINUED) 


Wheel 

Slip 

X 

Turn 

Angle 

a 

Degree 

Coefficients 

A 

B 

C 

20.0 

0 

0.640 

10.8 

3.81 

10 

0.556 

9.30 

3.5 

15 

0.365 

5.02 

4.0 

20 

0.221 

4.10 

7.3 

Equation  4.8  in  the  above  e3q>anded  version  was  used  to  formulate 
the  long-short  dashed  lines  Illustrated  In  Figures  4.14a  through  4.19a. 
Figures  4.14b  through  4.19b  illustrate  the  torque  coefficient  versus 
slip  at  various  sand  mobility  nuidjers  and  a range  of  wheel  turn  angles. 
Several  occurances  seem  apparent.  First  over  the  slip  range  tested, 
M/Wr  Increases  gradually  and  somewhat  llnearlly  as  slip  Increases  for 
a given  Ng  and  a . For  a given  wheel  slip  and  sand  mobility  number, 
M/Wr  decreases  as  the  wheel  turn  angle  Increases.  On  a percentage 
or  proportionate  bases  this  decrease  in  M/Wr  Is  not  as  large  as  that 

d 

noted  for  the  P/W  performance  parameter.  The  most  discerning  feature 
of  the  M/Wr  versus  wheel  slip  s plots  is  that  the  base  line  (i.e. 
a » 10)  predicted  by  Equation  2.10  frequently  He  below  the  plotted 
data  for  a turn  angle  of  five  degrees  and  on  occasion  plot  below  the 
data  points  corresponding  to  a 10  degree  turn  angle.  Examination  of 
the  basic  data  used  by  Tumage  and  Melzer  in  developing  M/Wr  versus 

A 

Ng  and  reported  herein  with  Figures  2,5  and  2,7  Indicates  that  for  Ng 
values  less  than  20  (the  majority  of  the  powered  turn  wheel  sand  tests 
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were  performed  with  less  than  20)  and  for  a given  wheel  slip,  the 

torque  coefficients  from  Individual  teats  varied  approximately  + 30  per- 
cent from  the  smooth  curve  drawn  through  the  data  points.  Comparisons 
of  the  data  points  of  Figure  A. 14b  through  4.19b  with  the  respective 
base  line  predictions  of  Equation  4.10  Indicates  that  more  than  80  per- 
cent of  the  M/Wr  test  data  Is  within  + 20  percent  of  that  value  pre- 
dieted  by  Equation  4.10  at  respective  values  of  wheel  slip.  Since  the 
difference  between  the  test  data  and  the  predicted  value  Is  less  than 
the  variation  of  the  basic  data  from  Equation  2.10  was  developed,  It 
does  not  seem  appropriate  to  refine  the  data  further  and  In  effect 
wheel  turn  angle  Is  considered  Insignificant  in  computing  torque  Inputs 
to  a powered  wheel  operating  In  loose  sand. 

Side  Forces  Developed  In  Sand 

Figure  4.28  illustrates  plots  of  the  side  force  coefficient  S/W 
versus  the  pull  coefficient  P/W  ; Individual  plots  contain  related 
data  having  about  equal  sand  mobility  numbers  Ng  . The  linear  rela- 
tions fitted  to  the  plotted  data  points  verify  the  finding  established 
for  the  clay  data;  namely,  that  for  positive  wheel  slip  values  up  to 
20  percent,  S/W  versus  P/W  for  a given  o and  Ng  can  be  expressed 
linearily.  Further,  the  slope  of  the  line  is  constant  for  a given  angle 
and  displaced  vertically  upward  on  the  plot  as  Ng  Increases.  The 
numbers  beside  the  plotted  points  Indicate  wheel  slip  & ; and,  as  with 
the  clay  data,  shows  that  as  the  Input  power  (as  expressed  and  implied 
by  wheel  slip)  to  the  wheel  Is  Increased,  the  pull  becomes  greater  and 
the  magnitude  of  the  side  force  decreases. 

Slopes  of  the  lines  drawn  through  the  data  points  in  Figure  4.28 
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Mobility  Nximber 


are  plotted  against  the  respective  turn  angle  a in  Figure  4.29.  A 
straight  line  was  passed  through  the  origin  and  the  plotted  points  hav- 
ing a slope  of  1.65  (note  that  this  value  is  only  slightly  below  the 
1.72  reported  on  Figure  4.6  from  similar  clay  data). 

The  values  of  side  force  coefficients  developed  at  zero  pull  were 
selected  from  Figure  4.28  and  plotted  in  Figure  4,30.  Figure  4.30  indi- 
cated that  thfe  side  force  coefficient,  when  the  pull  is  zero  and  the 
sand  mobility  number  is  about  the  same,  increases  with  Increasing  values 
of  turn  angle  a . Further,  if  the  turn  angle  is  held  constant,  S/W 
at  zero  pull  increases  as  the  sand  mobility  number  increases.  A second 
degree  polynominal  curve  was  fitted  to  data  points  having  Ng  of  8.3 
to  9.4  and  a second  polynominal  curve  for  those  tests  having  a Ng  of 
26.3.  These  curves  imply  that  Ng  is  of  secondary  importance  to  the 
development  of  side  forces  for  powered  wheels  operating  in  loose  sand. 

If  the  coefficients  of  the  two  polynominals  are  assumed  to  vary  linear- 
ily  for  intermediate  values  of  the  sand  mobility  number  then  a general 
equation  can  be  written  to  relate  the  side  force  coefficient  at  zero 
pull 

I = (2.3  + 0.03  N ) a - (2.4  + 0.065  N ) (4.9) 

= 0 

The  side  force  coefficient  at  zero  pull  is  the  ordinate  Intercept 
of  the  S/W  versus  P/W  relation  and  the  slope  1.65  a expresses  the 
slope  of  the  S/W  versus  P/W  relation.  Hence  for  a given  turn  angle 
and  sand  mobility  number  the  two  performance  parameters  are  related  as 

(«r ' (4 . 0 ■ w 
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SIDE  FORCE  COEFFICIENT  AT  P/W 


Figure  1|.30. 


Side  Force  Coefficient  at  P/W  = 0 Versus  Turn  Angle, 
a , for  Sand. 
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Figures  4.31  through  4.36  contain  the  plotted  test  values  of  side 
force  coefficient  versus  percent  wheel  slip.  Again  the  tests  are 
grouped  on  each  plot  having  similar  test  conditions  (i.e.  sand  mobility 
number  approximately  the  same) . Superimposed  are  broken  lines  deter- 
mined from  Equation  4.10  along  with  solid  lines  visually  fitted  to  the 
vUta.  For  the  sand  tests  conducted  with  turn  angles  other  than  zero, 
Equation  4.10  predicted  the  side  force  within  20  percent  or  less  of  the 
test  value  in  40  and  our  46  tests  or  83  percent. 

Slnkage  in  Sand 

Sinkage  coefficient  (z/d)  as  a function  of  the  sand  mobility  number 
and  wheel  turn  angle  is  Illustrated  in  Figure  4.37.  As  expected,  z/d 
decreases  for  a given  turn  angle  with  increasing  Ng  . The  broken  lines 
on  each  plot  of  Figure  4.37  represents  results  from  previous  test  pro- 
grams where  the  turn  angle  was  zero  and  the  wheel  slip  about  20  percent 
(Turnage,  1972).  Data  obtained  with  wheel  turn  angles  of  10,  15,  and 
20  degrees  indicate  that  for  a given  sand  mobility  number,  sinkage 
increases  with  increasing  turn  angle.  This  is  readily  apparent  from  the 
summary  plot  provided  in  Figure  4.37.  The  number  beside  each  datum 
point  is  the  wheel  slip  in  percent,  and,  as  with  the  clay  results,  sink- 
age  is  not  significantly  influenced  by  w)ieel  slips  up  to  20  percent. 

Summary 

For  a given  set  of  independent  variables  expressed  as  the  sand 
mobility  number  , the  wheel  slip,  and  the  wheel  turn  angle  forces 
acting  on  a pneumatic  tire  are  computed  as  follows: 

The  input  torque  coefficient  in  the  plane  of  the  wheel 
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Figure  I*. 33.  Influence  of  Turn  Angle  on  Side  Force  Coef f icient . I 

Average  Cone  Penetration  Resistance  Gradient  f 

G = 3.2  , Tire  Deflection  6/h  = 0.25  , Design  j 

Wheel  Load  W = 2,000  N . Desired  Sand  1 

Mobility  Number  N„  = 9.^  . ! 
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SIDE  FORCE  COEFFICIENT  S/W 


WHEEL  SUP  S,  PERCENT 

Figure  li.35.  Influence  of  T^rrn  Angle  on  Side  Force  Coefficient. 

Average  Cone  Penetration  Resistance  Gradient 

G = 2.0  , Tire  Deflection  6/h  = 0.35  j Design 

Wheel  Load  V;  = 1,000  N . Desired  Sand  Motility 

Numter  N = l6.4  . 
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SIDE  FORCE  COEFFICIENT  S/W 


Figure  I1.36.  Infliience  of  Turn  Angle  on  Side  Force  Coefficient. 

Average  Cone  Penetration  Resistance  Gradient 
G = 3.2  , Tire  Deflection  6/h  = 0.35  > Design 
Wheel  Load  W = 1,000  N . Desired  Sand  Motility 
Number  N = 26.3  . 


Wr  I Wr^ 


Ng  - C + B 


where 


A ■ 0.66 


B - 4.71  + 1.72/s 


C « -10 


The  pull  coefficient  In  the  plane  of  the  wheel 


„\o.  A B 

P.\  A a a 

Wy  a " -‘c  + B 

. / S a a 


where 


- Ap  - 1.42  a 


- Bp  - 16 


Ap  - 0.69  - 


Bp  - 10.8 


3I./3 


(0.6  + 3.1  s) 


and  finally  the  side  force  coefficient  perpendicular  to  the  plane  of  the 


wheel 
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CilAPTER  V 


RELATIONS  APPLIED  TO  A SIMPLIFIED  MODEL 

Control  and  Stability  of  a Four  Wheeled 
Vehicle  in  a Flat  Turn 


Introduction 

The  directional  stability  and  control  of  a four  wheeled  vehicle 
operating  on  soft  soils  is  studied  by  means  of  a simplified  theoretical 
analysis  which  takes  into  account  the  variation  of  the  cornering  per- 
formance of  pneumatic  tires. 

The  problem  of  directional  stability  and  control  in  a flat  turn  is 
formulated  with  steady-state  dynamical  equations  of  motion  having  two 
degree  freedom,  namely,  vehicle  yaw  and  vehicle  side  slip.  Rolling 
motions  of  the  spring  .mass  are  included  by  being  superimposed  on  the 
steady-state  analysis  in  order  to  enable  calculation  of  the  change  in 
vertical  loading  on  the  tires  resulting  from  vehicle  roll. 

Equations  of  Motion 

The  following  assumptions  were  made:  a)  parallel  tracks  for  both 

wheels;  b)  side  force  resulting  from  wheel  camber  insignificant;  c)  no 
lateral  tire  deformation;  d)  the  steering  wheel  is  held  fixed  at  a 
particular  setting,  and  the  steering  linkage  is  rigid.  Vehicle  roll  was 
not  included  as  a motion  coordinate,  but  weight  transfer  resulting  from 
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the  roll  about  the  loagltudinal  axis  of  a turning  vehicle  will  be 
discussed  la  the  next  section. 

Let  (x,y)  be  a set  of  Cartesian  axes  whose  origin  In  fixed  at  the 
vehicle’s  center  of  gravity  (c.g.).  The  x-axls  Is  along  the  longitu- 
dinal axis  of  the  vehicle  with  a positive  sense  toward  forward  motion. 

The  y-axls  runs  laterally  from  the  vehicle  c.g.  with  a positive  direc- 
tion toward  the  center  of  the  radius  of  curvature  along  which  the 
vehicle  travels.  The  positive  sense  of  the  axes  along  with  symbols 
of  pertinent  variables  are  shown  In  Figure  5.1. 

The  steady  state  response  of  the  vehicle  is  the  final  condition 
of  motion  of  the  vehicle  which  occurs  at  some  finite  time  after  the 
start  of  maneuver.  Here  It  Is  supposed  that  a specific  steer  angle  6 
is  applied  to  the  steered  wheels  and  held.  For  this  assessment  it  is 
sufficient  to  assume  that  the  angle  y between  the  constantly  applied 
drawbar  pull  Z and  the  direction  of  the  longitudinal  axis  of  the 
vehicle  is  constant. 

Suppose  that  the  c.g.  of  the  vehicle  is  moving  with  a constant 
forward  speed  V as  shown  in  Figure  5.1.  If  the  side  slip  angle  of 
the  vehicle  is  g , with  positive  convention  as  shown  in  Figure  5.1, 
then  for  small  values  of  B the  component  of  velocity  along  the  x-axis 
is  V*co8  B or  approximately  V ; in  the  y-axls  direction,  the  side 
slip  velocity  is  V*sin  3 or  approximately  V*3  . 

The  equilibrium  of  moments  in  terms  of  the  moving  Cx>y)  axes  read: 

I \p  » (a«cos  6 + c*sin  6)  + (a. cos  6 - d*sin  6)  - b (82^^ 

+T11  (c'cos  6 - a* sin  6)  +T^2  (d'cos  6 + a* sin  6)  + d T^2  ~ c T^^^ 

- Z L (5.1) 

y b 

With  the  instantaneous  yawing  velocity  notated  as  ip  and  positive 
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Forces  on  a Vehicle  Describing  a Circular 
llustration  of  the  Side  Slip  Angle. 
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as  skovm  in  Figure  5.1.  For  realistic  maneuveriag  of  vehicles  in  off- 
road terrains  steering  angle  6 is  sufficiently  small  that  cos  6»1 
and  sin  6 fed  . Equation  5.1  is  then  rewritten 

I (a  + cd)  + Ca  - dd)  - b CS21  + S22) 


(c  - a«  (d  + aS) 

(5.2) 

Tlie  equilibrium  of  forces  in  the  traverse  (y  axis)  and  longitudi- 

Vial  ( 

X axis)  directions  with  respect  to  the  vehicle  yield: 

M*V  ($  + ijj)  = cos  d + Sj^2‘aos  d + + 8^2 

+(T,  *sin  d + T._*sin  d)  + Z 
11  12  ' y 

(5.3) 

M-V  (B  + « - + S21  + S22  + (T^-6  + Ijj-S)  + 

(5.4) 

and 

’’^21  ’^22  + '■'12  «)  - 2^  * 0 

(5.5) 

"21  ^22  ’^12>  - "x  - » 

Also 

note  from  Figure  5.1  that 

Z = Z cos  Y 

X ' 

(5.6) 

and 

Z ■ Z sin  Y 

y ' 

(5.7) 

Those  terns  with  double  signs  the  upper  sign  signifies  a negative 
driving  force  (rear  wheel  drive)  at  the  front  wheels  and  the  lower  sign 

for  a positive  driving  force  (four  wheel  drive) . 


".Vi  • 
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In  a steady-state  turn  the  yawing  and  side  slip  accelerations 
become  zero  and  the  yawing  velocity  is  the  rate  of  turning  or 

e = 0 
^ = 0 

V cos 

♦ E E 

for  small  values  of  6 . 

Effective  Side  Slip  Angle 

The  forward  velocity  in  the  x-direction  at  the  center  of  each 
wheel  hub  is  equal  to 

Vu  - V cos  6 - c*  y 5) 

V^2  ~ V cos  E - di{i 

Recalling  the  assumption  that  g is  small  and  V cos  E«V  . Also  for 
steady  state  analysis  V = where  R is  the  radius  of  curvature  and 

R » c or  d . Compared  to  V , c*ij)  or  d*tj)  are  small  and  neglected. 

On  this  basis,  the  forward  velocity  in  the  x direction  at  each  wheel  is 
simply  V to  a sufficiently  close  approximation. 

The  lateral  or  side  slip  velocity  of  the  right  side  tires  arc 
schematically  represented  in  Figure  5.1,  and  equal  to  gV  + a'l'  for  the 
steered  tires  and  EV  - bi}»  for  the  rear  tires.  These  velocity  vectors 
form  an  angle  with  the  longitudinal  axis  of  the  front  and  rear  wheels 

3 + ■^  and  3 - ^ , respectively. 

V 

The  positive  sense  for  the  steer  angle  6 is  shown  in  Figure  5.1 
and  from  inspection  the  angle  between  the  steered  wheel  hub  and  the 
direction  of  travel,  designated  as  the  effective  side  slip  angle  a„  , is 

r 
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(5.9) 


a 


F 


6 - g 


ail) 
~ V 


The  side  forces  given,  in  the  equations  of  motion  (5.2  and  5.4)  have  their 
lines  of  action  normal  to  the  wheel  hub;  however  the  magnitude  of  these 
forces  depend  primarily  upon  the  angle  formed  by  the  side  slip  velocity 
vector  and  the  plane  of  the  wheel  hub.  Similarly  the  effective  side 
slip  angle  for  the  rear  wheels  can  be  similarly  shown  to  be 


a 


R 


(5.10) 


Vertical  Wheel  Loading 


To  maintain  the  vehicle  in  an  ideal,  steady-'State,  flat  turn,  the 
spring  mass  must  roll  toward  the  outside  of  the  curve  by  a certain 
amount  and  a subsequent  welglit  transfer  will  occur. 

To  compute  the  roll  angle  e for  a four  wheel  vehicle  describing 
a circular  curve  let  the  spring  mass  be  connected  to  wheels  with  springs 
having  elastic  constants  of  kp  and  corresponding  to  front  and 
rear  springs  respectively. 

Figure  5.2a  shows  the  representative  forces  and  sign  convention  of 
a four  wheeled  vehicle  undergoing  a steady-state  turn  of  radius  R at 
a constant  velocity  V . Figure  5.2b  is  a dynamical  equivalent  repre- 
sentation of  the  vehicle  with  the  roll  angle  having  been  developed  and 
the  wheels  and  suspension  replaced  with  idealized  springs. 

Summary  moment  about  the  e.g.,  about  which  roll  will  occur  and 

referring  to  Figure  5.2b  for  definition  of  terms; 

•• 

I>o  - Ic  = 0 for  steady  state  conditions 

2 

c^i  (kp  + V + d (kp  + kj^)  - h-2|-  (5.11) 
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vtiere  K is  the  vertical  distance  oi  the  vehicle  c.g.  above  the  ground. 
Note  that  roll  ang'  a e is  equal  approximately  (again  using  small  angle 
descriptions) 


(5.12) 


or 


» ce 


^2  ■ de 


(5.13) 

(5.14) 


Substituting  equation  5.13  and  5.14  into  equation  5.11  yields 

„2 


e = 


W 


c^  + d^ 


(5.15) 


Returning  to  Figure  5.1  the  vertical  weight  acting  through  the  c.g. 


is 


”cg  “ «11  ^ ”l2  + «21  + «22 


(5.16) 

The  static  weight  of  each  wheel  is  proportional  to  the  location  from  the 


c.g. 


W,,  = -W 


Wo 


^W. 


f 

12 

(a) 

“l2  ■ t “22 

(c) 

22 

(b) 

“22  ■ f “ll 

(d) 

(5.17) 


Repeated  substitution  of  equation  5.17  into  equation  5.16  yields  the 
static  weight  at  each  wheel 

bd 


% * (a  + b)  (c  + d)  '^cg 

12  (a  + b)  (c  + d)  '^cg 


(5.18) 


W 
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^21  (a  + b)Cc  + d)  '^cg 


^^22  ~ (a  + b)  (c  + d)  '^cg 


Next  the  vertical  force  at  each  wheel  during  steady-state  turning 
of  the  vehicle  is  computed  by  issuming  the  forces  to  be  composed  of  the 
static  component  weight  acting  when  the  vehicle  is  at  rest  on  a level 
horizontal  plane  plus  the  component  arising  as  a result  of  dynamic  body 
roll.  The  weight  at  each  wheel  will  then  be 


W..  = W. . +k.  ex 
ij  ij  - j 


(5.19) 


where  x is  the  left  or  right  vehicle  dimension.  The  upper  sign 
designates  the  outside  wheels  and  the  lower  sign  for  the  inside  wheels. 
As  previously  subscript  i refers  to  front  (i  * 1)  and  rear  (i  = 2) 
axle  and  j designates  whether  the  wheel  is  nearest  the  inside  of 
curvature  ( j * 1) . 


W = Siit 

11  (a  + b)  (c  + d) 


W s=  — 

12  (a  + b)(c  + d) 


21  (a  + b) (c  + d) 


22  (a  + b)  (c  + d) 


W - kp  ce 

(5.19a) 

W + k„  de 

(5.19b) 

W - kp  ce 

(5.19c) 

W + kp  de 

K 

(5.19d) 

Tire  Deflection  Dependency  Upon 
Vflieel  Load  and  Tire  Inflation  Pressure 


For  a four  wheel  vehicle  undergoing  a steady-state  turn  the  load 
for  each  wheel  and  the  corresponding  tire  deflection  probably  differs. 
The  clay  and  sand  numerics  introduced  and  discussed  in  Chapter  II  and 


' * ' ^ . 4 V 


utilized  in  Chapter  IV  vary  inversely  with  changes  in  load  W , and 
either  directly  or  proportionately  as  the  deflection  ratio  6/h 
fluctuates.  If  ■'  change  in  these  variables  occur  a procedure  is 
required  to  ensure  that  the  numeric  property  reflect  existing  para- 
meters. Numerous  tires  have  undergone  static  dynometer  testing  at 
the  Mobility  Testing  Facility  at  the  WES  in  establishing  interrelated 
data  on  wheel  load,  tire  deflection,  and  inflation  pressure.  Figures 
5.3  and  5.4  illustrate  the  seemingly  linear  dependency  of  the  ratio 
of  wheel  load  W to  tire  deflection  ratio  6/h  on  inflation  pressTire 
IP  . Figure  5.3  contains  W/(6/h)  versus  IP  data  for  the  6.00-9, 

4-PR  tire  used  in  the  test  program  and  discussed  in  Chapter  IV. 

Figure  5.4  illustrates  W/(6/h)  versus  IP  data  for  four  times  common 
to  U.  S.  Army  tactical  wheeled  vehicles. 

A linear  equation  based  on  least-squares-fit  has  been  developed 
for  each  set  of  W/(6/h)  versus  IP  data.  Having  tire  deflection 
expressed  as  equations  continuous  across  possible  range  of  tire  infla- 
tion pressures  is  inqjortant  because  as  roll  and  subsequent  weight 
transfer  occur  a change  of  the  tire  deflection  results.  Hence  if 
inflation  pressure  is  considered  constant^  and  the  load  per  wheel  is 
computed  from  the  roll  equation  of  the  last  section  then  the  resulting 
deflection  can  be  computed  from  the  corresponding  equation  provided 
on  Figures  5.3  and  5.4. 


Constant  inflation  pressure  of  military  tires  is  a meaningful 
assumption  since  the  pressures  for  either  on— road  or  off-road  operations 
is  .‘generally  designated  on  the  vehicle  and  strictly  enforced. 
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Figure  5.3.  Dependency  of  Load/Deflection  Rata.o  on  Tire 
Pressiire. 
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800 
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INFLATION  PRESSURE  IP,  KN/M^ 

Figure  5.4.  Dependency  of  Load/Def lection  Ratio  on  Tire  Pressure 
and  Size.  ^ 


Turned  Tlre/Soll  Interaction  Application 


PlKltal  Model 

The  dyntvnlcal  degrees  of  freedom  chosen  for  analysis  consist  of 
vehicle  yaw  and  vehicle  side  slip.  Vehicle  roll  is  not  considered  as  a 
dynamical  coordinate  but  as  a quasi-static  factor  dependent  upon  the 
yaw  velocity.  With  the  axes  fixed  in  the  body  of  the.  vehicle  the  two 
degree  of  freedom  model  may  be  expressed  by 

MV  (P  + i|j)  = F 

y 

and 

= M 
^ z 

which  are  notational  forms  of  Equations  5.5  and  5.2,  respectively. 
Solving  the  above  set  of  equations  require  that  tractive  and  side 
forces  be  determined,  which,  in  turn,  are  dependent  upon  vehicle  para- 
meters and  soil  conditions  in  addition  to  nonlinear  functions  of  wheel 
side  slip.  A alogrithm  developed  by  Brown  (1967)  was  adopted  to  solve 
the  two  simultaneous  nonlinear  equations  having  unknowns  of  vehicle 
yaw  ip  and  vehicle  slip  6 . 

Brown’s  alogrithm  is  a modification  of  Newton’s  method  for  solving 
simultaneous  nonlinear  equations,  requiring  no  derivative  evaluations. 

As  with  most  iterative  numerical  techniques,  the  closer  the  initial 
guess  is  to  the  true  solution  the  fewer  iterations  that  will  be  neces- 
sary for  a closed  solution;  also  the  accuracy  or  even  the  ability  of 
the  alogrithm  to  converg  toward  the  solution  often  depends  upon  the 
initial  guess  of  the  variables.  To  aid  in  selecting  a close  apprpxima- 
tion  of  and  3 for  initial  input  values  to  the  Brown’s  alogrithm,  a 
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method  was  adopted  that  is  used  for  finding  complex  zeros  of  functions 
(Hamming,  1973). 

The  previous  equations  can  be  rewritten  as 

MV  e * iy  - MVip 
and 

In  a steady  state  turn  the  yawing  and  side  slip  accelerations,  tj)  and 
B , respectively,  are  zero  making  the  right  side  of  the  above  equations 
equal  zero  when  compatible  values  of  6 and  ij)  are  known.  Finding 
real  zeros  of  the  above  equations  can  be  very  easy  and  yet  very  robust. 

( 

{ 

The  approach  is  to  assume  a set  of  values  for  3 and  tj/  and  solve  • 

1 

Equations  5.2  and  5.5  arranged  to  the  above  form.  Yaw  velocity  and 

vehicle  side  slip  values  are  varied  about  a m x n matrix  in  which  m ' 

represent  the  number  of  tj)  values  m\d?.r  consideration  and  n repre-  ; 


sents  the  number  of  3 values.  Generally  ip  is  varied  from  0 to  0.6 
radians  per  second  and  3 is  varied  from  minus  to  plus  1.5  times  the  . 

wheel  steer  angle  6 at  each  location  of  the  matrix,  values  of  3 and  i 

6 are  used  in  Equations  5.2  and  5.5  and  a value  of  1,  2,  3,  or  4 is  ; 

recorded  at  each  point  in  the  matrix  depending  upon  the  sign  convention  \ 

of  I'{i  and  MV  3 as  illustrated  in  Figure  5.5.  Obviously  when  the  - 

i 

values  arc  printed  for  the  m x n.  matrix,  where  the  four  quadrants  meet  \ 

I 

1 

at  a point  iu  the  general  region  of  zero.  The  solution  could  be  made  J 

more  refined  by  enlarging  the  zero  region  with  a finer  and  finer  grid  j; 

spacing.  However,  the  system  requires  too  much  interaction  with  the  ‘ 

i 

operator  and  after  an  approximate  determination  of  and  3 are  ; 

known  then  the  Brown's  alogritbjn  can  be  effectively  used  to  refine  the  i 
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solution.  Generally  a matrix  size  of  60  x 60  is  used  since  this  size 
is  readily  adaptable  to  a computer  time  sharing  teletype  terminal. 


+ MV  3 


Figure  5.5.  Scheme  for  Finding  Real  Zeros. 


Tlie  flovr  diagram  for  determining  approximate  values  of  tjj  and  3 
is  shown  in  Figure  5.6.  The  input  data  are  read  in  by  data  statements. 
Required  input  data,  designation  of  variable  names,  and  engineering 
units  are  listed  in  Table  V. 


TABLE  V 

LISTING  OF  INPUT  DATA  FOR  COMPUTER  MODEL 


Name 


Description 


Dimension 


Vehicle 

WIG  Gross  vehicle  weight  Newtons 

ADMIN  Distance  between  front  axial  and  vehicle  CG  Metre 

BDIM  Distance  between  rear  axial  and  vehicle  CG  Metre 

CDIM  Horizontal  distance  between  .<nside  wheels  and  Metre 

vehicle  CG 

DDIM  Horizuital  distance  between  outside  wheels  and  Metre 

vehicle  CG 

HDIM  Vertical  distance  from  ground  surface  to  Metre 

vehicle  CG 
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TABLE  V (ContJjiued) 


Name 

Description 

Dimension 

Vehicle 

LDIM 

Distance  point  of  application  of  drawbar  pull 

Metre 

to  CG 

GAMMA 

Angle  between  drawbar  pull  and  longitudinal 

Degrees 

axis  of  vehicle  2 

CTR 

Ratio  of  centrifugal  acceleration  (V  /R)  to 

— 

gravity  acceleration  (G) 

Newtons /Metre 

SPKE 

Spring  constant  for  front  axle 

SPKR 

Spring  constant  for  rear  axle 

Newtons /Metre 

Tire 

SLIP 

Powered  wheel  slip 

— 

TIRED 

Cross  section  width  of  unloaded-inflated  tire 

Metre 

TIRED 

Diameter  of  unloaded-inflated  tire 

Metre 

DEFI 

Tire  deflection 

— 

ITIR 

Tire  code  - from  MSD  tire  book 

— 

TRIP 

Tire  Inflation  pressure 

kPa 

Soil 

KSOIL 

0 ■ clay;  1 » sand 

— 

CPR 

Cone  penetration  resistance,  clay 

kPa 

CPR 

Cone  penetration  resistance,  sand 

MPa/m 

Parametric  Variation  for  Specific  Vehicle 

KDRIVE 

+1  ■ rear  wheel  drive;  -1  * four  wheel  drive 

MM 

VEL 

Vehicle  speed 

Metre/Second 

DELTA 

Steering  angle  of  steered  wheels 

Degrees 

lEND 

0 » more  data;  1 * last  data 

Approximate  values  of  if/  and  3 are  then  placed  as  input  to  a 
driver  computer  program  which  has  a subroutine  utilizing  Brown's 
alogrithm  and  another  subroutine  to  solve  Equations  5.2  and  5.5.  The 
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flow  diagram  is  illustrated  iu  Figure  5.7  aad  the  complete  computer 
program  listed  in  Appendix  C. 

The  computation  of  thrust  and  side  forces  is  the  same  irrespec- 

2 

tive  of  whether  the  program  for  finding  real  zeros  or  the  general 
turning  program  utilizing  the  Brown's  alogritlvm  is  used.  Significant 
computations  begin  by  determining  front  and  rear  wheel  slip  angles 
from  Equations  5.9  and  5.10,  respectively,  from  input  values  of  forward 
velocity,  fore  and  aft  displacement  of  the  axles  from  the  vehicles 
center  of  gravity,  and  assumed  values  of  vehicle  yaw  and  side  slip 
angle.  Next  the  program  computes  the  vehicle  roll  about  the  x-axis 
and  respective  wheel  loads  per  Equations  5.19a  through  5.19d.  A mobil- 
ity number  is  computed  for  each  wheel  using  Equation  2.2  if  the  surface 
soil  is  sand  and  Equation  2.3  for  a soft  clay  medium.  If  either  of  the 
clay  mobility  numbers  is  less  than  7 or  else  either  of  the  sand  mobil- 
ity numbers  is  less  than  2;  then  immobility  of  the  vehicle  is  said  to 
exist.  Examination  of  dependent  performance  parameters  illustrated  in 
Figure  2.4  for  clay  and  Figure  2.5  for  sand  indicates  that  for  N^=  2.5 

(approximately  equivalent  to  N*  = 7)  and  N = 2 that  sinkage  and 

c s 

resistance  to  pull  and  tow  become  excessive  whereas  forward  movement 
would  cease.  The  program  branches  depending  upon  whether  the  soil  is 
sand  or  clay  and  upon  whether  rear  drive  or  four  wheel  drive  is  employed 
for  computing  tractive  and  side  forces  acting  at  each  wheel.  Finally 


A listing  of  the  computer  program  for  finding  real  zeros  is  not 
given  since,  except  for  the  DO  loops  used  to  step  through  the  m x n 
matrix  and  output  print  statements  signifying  the  sign  of  and  MV  S 
both  amply  illustrated  in  Figure  5.6,  it  is  identical  to  the  sub-  ^ 
routine  TUEN  used  in  the  computer  program  listing  in  Appendix  C. 


- 

t . o , 
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these  forces  and  their  respective  moments  are  summed  according  to 
Equations  5.2  and  5.5,  and  if  the  solution  of  the  system  has  sufficient- 
ly converged,  output  of  pertinent  forces,  velocities,  and  turning 
configuration  are  displayed  to  the  user. 

All  of  the  horizontal  components  of  force  acting  on  a pcjwered 
wheel  and  in  the  direction  of  travel  can  be  summed,  and  their  sum  de- 
fined as  thrust  T . Tlie  sum  of  the  horizontal  forces  acting  parallel 
with  the  wheel  hub  and  opposite  to  the  direction  of  travel  can  be 
defined  as  rolling  resistance.  From  equilibrium  conditions,  thrust  is 
equal  to  rolling  resistance  plus  pull.  Because  both  thrust  and  rolling 
resistance  occur  in  the  same  physical  area,  they  are  difficult  to 
separate;  however  their  difference,  pull,  can  be  measured.  The  input 
torque  M to  the  powered  wheel  can  also  be  measured  and  it  represents 
the  thrust  times  the  lever  arm  (the  distance  below  the  axle  where  the 
horizontal  thrust  acts).  Also  the  difference  between  the  torque  and  pull 
coefficient  at  any  positive  slip  (see  Figure  2.2)  represents  resistance 
losses.  Besides  rolling  resistance  there  are  also  mechanical  (mainly 
frictional)  losses  of  the  wheel  and  forces  expended  in  deforming  the 
pneumatic  tire.  However,  the  tires  used  to  develop  all  relations  pre- 
sented thus  far  have  side  walls  which  flex  rather  easily;  hence  that 
force  required  for  deforming  the  pneumatic  tire  is  small  for  relations 
developed  herein.  Also  the  laboratory  wheel  bearing  has  been  well 
designed  and  maintained  to  reduce  friction  to  a minimum.  Therefore  for 
purposes  of  this  study  the  torque  coefficient  was  taken  as  the  thrust 
coefficient. 
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Sample  Calculations 


To  illustrate  the  nature  of  the  results  which  follow  from  the 
analysis,  numerical  calculations  have  been  carried  out  and  the  changes 
of  predicted  forces  and  yawing  vehicles  described  wiien  steady-state 
flat  turns  are  traversed  over  a sand  of  uniform  consistency  at  various 
speeds  and  curvatures.  In  the  following  numerical  example  the  steady- 
state  condition  for  a rear-wheel-drive  and  a four-wheel-drive  vehicle 
was  calculated  corresponding  to  the  previously  derived  expressions  and 
limiting  assumptions.  The  vehicle  data  presented  in  Table  VI  cor- 
responds closely  to  a M-151,  1/4- ton  capacity  military  jeep. 


TABLE  VI 

VEHICLE  DATA  USED  IN  NUMERICAL  EXAMPLE 


Item 

Symbol 

Units 

Gross  vehicle  weight 

W 

14,240  N 

Drawbar  pull 

Z 

0 N 

Tire  size 

— 

7.00-16,  6 PR 

Tire  pressure 

IP 

103.5  kN/M2 

Wheel  slip 

. S 

20% 

Front  axle  spring  constant 

25,228  N/M 

Rear  axle  spring  constant 
Vehicle  dimensions  (see  Figure  5.1) 

25,228  N/M 

a 

— 

1,143  M 

b 

1.016  M 

c 

— 

0.890  M 

d 

0.890  M 

h. 

— 

0.629  M 

Z 

— 

1.758  M 

V 

0® 

Steer  angle 

6 

7° 

Sand  penetration  resistance 

G 

2.5  MPa/m 
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The  steer  angle  was  kept  constant  and  the  forward  speed  of  the 
vehicle  was  the  variable  parameter.  Figure  5.8  provides  a comparison 
of  path  curvature  1/R  = <|)/V  , yaw  velocity  r)»  , vehicle  side  slip 
angle  8 , and  vehicle  body  roll  e versus  speed.  The  advantage  of 
four-wheel-drive  over  rear-wheel  drive  when  pulling  on  a curve  is  seen 
to  become  more  apparent  with  increases  of  forward  velocity  from  the 
1/R  versus  V plot  of  Figure  5.8.  As  V approaclves  zero  R ap- 
proaches Ackermann  neutral  steer  which  can  be  closely  approximated  by 


a + b 


neutral 


= 17.7  m 


or  for  path  curvature,  1/R  = 0.057  m~^.  The  smaller  turning  radius  at 
equal  speeds  exhibited  by  the  example  vehicle  in  the  four-wheel  drive 
mode  is  also  reflected  in  higher  yaw  velocity  and  greater  vehicle  body 
roll  as  conq)ared  to  the  rear-wheel-drive  configuration.  The  vehicle 
side  slip  angle  begins  at  an  angle  equal  to  the  steer  angle  as  the 
velocity  approaches  zero  and  has  a positive  sense  the  same  as  that  of 
the  steered  front  wheels.  As  speed  increases  8 reduces  such  that  the 
velocity  vector  approaches  the  longitudinal  axis  of  the  vehicle  and 
with  further  increases  in  forward  speed  8 again  increases  in  magnitude 
but  in  the  opposite  angular  direction  from  that  of  the  steered  wheels. 
Over  rhe  range  of  speeds  investigated,  the  vehicle  side  slip  angle 
changes  considerably  faster  for  the  four-wheel-drive  vehicle. 

A pictoral  comparison  is  presented  in  Figure  5.9  of  the  equilib- 
rium configuration  tliat  are  established  in  a right  turn  for  the  rear- 
driven  and  all-wheel- drive  at  4-m/sec  and  steer  angle  of  7 degrees. 
Included  on  each  diagram  are  the  three  components  of  force  on  each  tire, 
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Figure  5*8.  Comparison  of  Motion  Inducted  Parameters  Versus  Velocity 
for  Rear  VRieel  and  Four  Wheel  Driven  Vehicle.  Steer 
Angle  6=7  Degrees  . 
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and  the  vehicle  responses  V , 'I*  , P > Ny  Ccentrifugal  acceleration 
force),  and  R . The  proper  angular  orientation  and  magnitudes  of 
the  horizontal  forces  end  velocity  vectors  are  shown.  Note  that  for 
the  rear-drive  vehicle  Oj,  > and  R > \£^^.ral’ 

all-wheel  driven  a^,  < and  R < would  indicate  by 

definition  that  at  V = 4 m/s  and  6=7  degrees  the  vehicle  is  in  an 
oversteer  condition  when  all-wheel— drive  is  utilized  and  in  an  under- 
steer condition  when  only  rear-drive  pervails.  Tliese  conditions  of 
oversteer  and  understeer  over  the  full  velocity  range  considered  are 
evident  from  the  1/R  versus  V relation  of  Figure  5.8  in  that  1/R 
increases  with  increases  of  V for  the  rear  driven  version. 


CHAPTER  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 
Conclusions 

Based  on  the  analysis  of  the  data  obtained  during  this  research 
the  following  conclusions  can  be  drawn  with  regard  to  the  performance  of 
single  powered  wheels  equipped  with  pneumatic  tires  and  operating  in  the 
turn  mode  on  soft  yielding  soils. 

1)  For  the  test  variables  considered  in  this  study,  the  WES 
systems  for  predicting  the  performance  of  powered  wheels  operating  at 
zero  turn  angle  on  saturated  clay  and  air-dried  sand  can  be  extended 
and  modified  to  treat  turn  angles  larger  than  zero  and  wheel  slips  other 
than  20  percent  such  that  the  principal  performance  parameters  of  pull 
and  torque  can  be  effectively  predicted. 

2)  Performance  expressed  in  terms  of  pull  coefficient,  side  force 
coefficient,  input  torque  coefficient,  and  sinkage  coefficient  was  found 
to  be  influenced  by  the  same  independent  variables  identified  in  pre- 
vious studies  plus  that  of  wheel  turn  angle.  From  the  test  results  the 
following  generalities  can  be  made  concerning  the  individual  relations 
of  the  performance  parameters; 

a)  At  specific  values  of  wheel  slip  and  mobility  number, 

increasing  the  wheel  turn  angle  reduces  the  pull  coefficient. 
For  a given  clay  mobility  number,  the  reduction  in  pull 
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coefficient  from  that  associated  with  a zero  turn  angle 
for  the  powered  turned  wheel  operating  with  a constzmt 
turn  angle  is  nearly  constant  over  the  wheel  slip  range 
tested;  however,  for  tests  performed  on  sand  the  amount  of 
decrease  of  the  pull  coefficient  lessens  with  increases  in 
wheel  slip. 

b)  It  was  showii  that  by  increasing  the  driving  force  at  con- 
stant wheel  turn  angle  and  mobility  number  the  side  force 
was  considerably  reduced.  For  a specified  wheel  turn  angle 
and  mobility  number  the  relation  between  side  force  coef- 
ficient and  pull  coefficient  can  be  approximated  with  a 
straight  line  over  the  wheel  slip  range  tested.  This 
observation  was  found  to  be  true  for  both  clay  and  sand 
tests. 

c)  If  the  independent  variables  expressed  in  terms  of  mobil- 
ity number  and  the  wheel  slip  is  kept  constant,,  the  side 
force  coefficient  increases  as  the  wheel  turn  angle  becomes 
greater.  Further,  if  the  wheel  turn  angle  and  slip  is  held 
constant,  the  side  force  coefficient  increases  with  in- 
creases in  the  mobility  number.  Finally,  if  the  wheel  turn 
angle  and  mobility  number  has  specific  values,  the  side 
force  coefficient  decreases  with  increases  of  wheel  slip. 

d)  Sinkage  coefficients  for  clay  and  sand  decreases  at  a power 
decay  rate  as  the  mobility  number  increases  for  a given 
wheel  turn  angle.  For  a given  sand  mobility  number,  the 
sinkage  coefficient  increases  with  increasing  wheel  turn 
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pngle.  Over  the  range  tested,  wheel  slip  does  not 
appreciably  effect  sinkage  in  either  clay  or  sand. 

Experimental  results  on  single  tires  running  at  various  turn 
angles  were  applied  via  equations  of  motion  to  the  analysis  of  four- 
wheel  vehicles  executing  a turn  on  flat,  level,  soft  soil.  The  steady- 
state  turning  behavior  of  four-wheeled  vehicles  at  different  values 
of  steer  angle  and  speed  was  studied  to  illustrate  how  the  interaction 
01  the  side  force  and  the  driving  force  affects  the  motion  of  the 
vehicle. 


Recommendations 

In  order  to  develop  a better  understanding  of  wheeled  vehicle 
agility  in  off-road  operations,  the  following  recommendations  for  fur- 
ther research  may  prove  useful; 

1)  Research  reported  in  this  thesis  is  limited  to  application 
of  wheeled  vehicles  in  soft  soil.  A study  would  be  beneficial  that 
investigated  the  turned  tire-soil  Interaction  when  operating  in  soils 
of  medium  to  stiff  consistency.,  Tliese  results  could  be  combined  with 
the  results  reported  herein  concerned  with  soft  soil  and  from  the 
numerous  published  results  of  turned  tires  operating  on  rigid  or  flexi- 
ble pavements  toward  developing  generalized  relations  for  tire-surface 
interaction. 

2)  The  example  used  in  Chapter  V to  illustrate  steady  state 
turning  behavior  indicates  that  effective  side  slip  at  the  wheels  are 
small  in  magnitude.  Hence,  laboratory  tests  should  be  extended  to 
examine  the  side  forces  developed  with  wheel  turn  angles  of  less  than 
5 degrees. 
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3)  Formulate  rationale  for  bridging  steered  tire-soil  inter- 
action relations  developed  from  laboratory  experiments  using  cohesion- 
less sands  and  a saturated  cohesive  clay  to  wheeled-vehicles  operating 
in  intermediate  soils  having  both  cohesive  and  frictional  properties. 

4)  Determine  the  influence  of  a thin,  slippery  surface  layer  of 
soil  on  the  performance  of  off-road  wheeled  vehicles  while  turning. 

5)  Conduct  laboratory  testing  program  to  determine  the  influence 
that  braking  has  on  the  development  of  side  forces  on  a turned  tire. 

6)  Incorporate  the  effects  on  vehicle  performance  produced  by 
terrain  side  slope  and  braking  into  the  computer  model  developed  during 
this  study  for  predicting  the  time  history  of  a wheeled  vehicle  travers- 
ing a prescribed  path.  This  path  following  model  should  be  developed 

to  exploit  as  much  of  the  existing  (straight  line  travel)  AMC  Mobility 
Model  philosophy,  modeling  technique,  and  terrain  representation  as 
possible. 
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appendix  a 


TABULATION  OF  EESULTS  FROM 
SINGLE  WHEEL  TIRE  TESTS 


(effect  of  trail  moment  neglected) 


appendix  b 


mobility  number  computed  from 
measured  test  values 


Clay  Tests 


The  dependent  paraaetexa  of  gull,  input  torque,  and  wheel  slip 
nveasured  in  each,  clay  test  can  be  substituted  into  Equations  2.4 
through  2.7  and  the  clay  mobility  number  computed  which  in  turn  can  be 
compared  with  the  clay  mobility  number  determined  by  test  conditions 
for  assessing  the  influence  of  wheel  turn  angle  on  performance  para- 
meters. Three  combinations  of  the  above  mentioned  measured  dependent 
parameters  can  be  used  for  computing  the  clay  mobility  number:  (1)  pull 

and  wheel  slip,  (2)  input  torque  and  pull,  and  (3)  input  torque  and 
wheel  slip. 


Pull  and  Wheel  Slip 


Beginning  with  Equation  2.5 

S 


21 


sp 

and  substituting  into  Equation  2.7 


P 1,  S 1, 

V - S~  ■ 2^°8 

sp 


+ 0.005 


21 


,x5/2 
c 


<K) 


+ 0.005 


and  finally  solving  for  N’ 


N' 

c 


21 


2/5 


I -1  1 

[log  \^W  y 


- 0.005 


(2.5) 


(B.l) 
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Input  Torque  and  Pull 
M 

Solving  for  from  Equation  2.4  yields 
*^a 


!k  » _ i (1  + b)  1/4 

Wr  Wr  W d'' 
a a 


(2.4) 


Equation  2.6  states 


M 

_££ 

Wr 

a 


2 


+ 0.007 


(2.6) 


Equating  the  right  side  of  the  two  previous  relations  and  solving  for 


yields 


N*  = 
c 


12 


- I (1  + Y'‘‘  - 0.007 
a 


1/2 


(B.2) 


Input  Torque  and  Wheel  Slip 


Smith  (1976)  found  that  the  torque  coefficient  and,  the  slip  at  the 
self-propelled  point  are  related  to  the  clay  mobility  number  by 


and 


M 

= 

Wr 

a 


+ 0.007 


S 

ap 


0.005 


(2.6) 


(2.5) 


Substituting  the  above  into  Equation  2,7, 


M - M 


Wr, 


SS. 


h 


C2.7) 
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creates  an  expression  having  terns  or  input  torque,  wheel  slip,  and  the 
clay  mobixity  number,  i.e. 


- log  S =«  + 0.014  - log  + 0.005  (B.3) 

a M N 

c c 

can  be  calculated  for  specific  values  of  input  torque  and  wheel  slip 
from  the  above  relation  by  a simple  interaction  algorithm  programmed 
into  a computer. 

Table  B.l  provides  a comparison  of  the  clay  mobility  number 
as  computed  from  Equations  B.l,  B.2,  and  B.3.  These  values  did  not 
differ  significantly  however,  as  computed  from  the  pull  and  wheel 
slip  test  values  were  used  in  further  developments  since  they  represent 
system  outputs  of  the  powered  wheel. 


TABLE  B.l 

CLAY  MOBILITY  NUMBER  COMPUTED 
FROM  PERFORMANCE  PARAMETERS 


Tes.t.Ngi, 


Angle 


Computed  Clay  Mobility  Number 
Pull Torque Torque 

& & S> 

Slip  , Pull  Slip  Avg. 


A-73-0U48-3 

0 

14.65 

A-73-0049-3 

0 

18.63 

A-73-0050-3 

0 

19.35 

A-73-0051-3 

10 

13.96 

A-73-0052-3 

10 

14.16 

A-73-0053-3 

10 

15.89 

A-73-0054-3 

20 

10.73 

A-73-0055-3 

.20 

9.67 

A- 73-0056-3 

20 

10.87 

A- 7 3-005 7-3 

5 

15.72 

14.52 

14.58 

14.58 

18.84 

13.99 

17.35 

16.66 

20.00 

15.59 

17.55 

17.50 

17.87 

11.92 

12.94 

12.94 

17.34 

11.74 

13.19 

13.03 

18.97 

12.86 

14.77 

14.51 

18.88 

29.94 

16.41 

19.03 

19.75 

12.67 

9.65 

10.67 

19.50 

26.10, 

14.57 

17.18 

19.88 

13.40 

14.56 

14.56 

18.61 

TABLE  B.l  (Continued) 


Computed  Clay  Mobility  Number 


Pull 

Angle 

& 

Test  No. 

a,  deg. 

Slip 

A-73-0058-3 

5 

11.29 

A-73-0059-3 

5 

10.46 

A-73-0060-3 

5 

9.29 

A-73-0061-3 

15 

12.08 

A-73-0062-3 

15 

8.43 

A- 7 3-006 3- 3 

15 

A-73-0064-3 

10 

12.49 

A-73-0065-3 

10 

12.90 

A-73-0066-3 

10 

li.52 

A- 7 3-006 7- 3 

10 

6.51 

A-73-0068-3 

10 

25.92 

A-73-0069-3 

5 

10.68 

A-73-0070-3 

0 

25.99 

Torque 

Torque 

& 

£t 

N' 

c 

Pull 

Slip 

Avg. 

10.92 

10.92 

12.59 

10.55 

12.63 

9.44 

9.36 

9.36 

13.00 

11.40 

11.74 

11.74 

11.96 

9.60 

8.93 

8.98 

12.06 

12.13 

8.20 

9.74 

12.93 

15.55 

13.38 

13.81 

15.44 

13.13 

13.02 

14.64 

11.51 

11.51 

14.23 

8.01 

3.34 

5.96 

7.54 

11.81 

22.78 

20.17 

30.32 

10.76 

10.72 

10.72 

12.60 

17.21 

23,52 

22.24 

18.97 

Sand  Tests 

Equations  2.9  and  2.10  were  used  to  calculate  the.  sand  mobi^-ity 
number  from  test  values  of  wheel  slip,  input  torque,  and  pull.  At  this 
tine  an  equation  derived  from  test  data  of  input  torque  and  wheel  slip 
at  the  self-propelled  has  not  been  developed;  hence,  there  is  one  less 
set  of  relations  for  the  sand  tests  from  which  to  compute  the  sand 
mobility  number  from  dependent  test  parameters.  Table  B.2  provides 
a comparison  of  the  sand  mobility  number  as  confuted  by  substituting 
combinations  of  dependent  performance  parameters  pull-wheel  slip  and 
torque-wheel  slip  into  Ec^uations  2.9  and  2.10,  respectively.  Also  shown 
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are  Ng  values  computed  from  Independent  test  conditions.  As  seen  in 
Table  B.2  the  computed  sand  mobility  number  was  highly  dependent  upon 
the  method  and  equation  selected.  Further,  when  compared  with  the 
respective  Ng  determined  by  test  conditions,  the  computed  values  did 
not  follow  any  well  defined  trend  as  values  occurred  almost  equally 
above  and  below  the  base  line  N_  value. 

O 


TABLE  B.2 

SAND  MOBILITY  NUMBER  COMPUTED 
FROM  PERFORMANCE  PARAMETERS 


Test  No. 


Wheel 
Turn 
Angle 
o.  deg. 


Computed  Sand 
Mobility  Number 
Pull  Torque 

£>  & 
Slip  Slip 


. A-73-0011-1 

0 

8.02 

15.27 

29.02 

i A-73-0012-1 

0 

7.85 

11.35 

20.44 

1 A-73-0013-1 

6 

7.62 

10.31 

17.66 

i A-75-0011-1 

5 

8.57 

•k 

* 

1 A-75-0012-1 

5 

8.10 

11.82 

18.56 

1 -A-75-0013-1 

5 

8.64 

10.44 

17.98 

1 A-73-0014-1 

5 

8.29 

8.78 

16.10 

i A-7A-0005-1 

10 

7.68 

12.62 

18.78 

A-74-0001-1 

10 

8.70 

9.16 

13.87 

A-74-0003-1 

10 

8.54 

8.18 

13.59 

A-74-0004-1 

10 

8.58 

9.24 

16.31 

! A-74-0002-1 

10 

8.59 

7.57 

11.45 

! A-73-0015-1 

10 

7.61 

7.41 

15.84 

A-74-0006-1 

20 

8.18 

4.85 

3.17 

A-74-0007~i 

20 

8.37 

4.60 

2.88 

A-74-0024-1 

15 

17.53 

31.22 

23.81 

A-75-0015-1 

20 

17.19 

n>  rs 

0.^0 

t /:  AT 

XU  .U  '/ 

A-75-0014-'l 

20 

16.08 

6.12 

11.72 

1 

) * Not  computed  because  wheel 

i : 

1 

slip  was  negative. 

- 
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TABLE  B.2  (Continued) 


Wheel 

Turn 


Computed  Sand 
Mobility  Number 


Angle 

N 

Sc 

S. 

Test  No. 

a,  deg. 

s. 

Slip 

Slip 

A-75-0016-1 

20 

17.43 

6.99 

15.05 

A-74-0023-1 

15 

5.04 

4.11 

5.11 

At  75-000 1-1 

10 

4.01 

19.55 

27.68 

A-75-0002-1 

10 

3.76 

4.64 

11.41 

A-75-0003-1 

10 

3.49 

3.58 

0.19 

A-74-0008-1 

20 

3.49 

4.61 

6.47 

A-74-0025-1 

10 

5.85 

6.24 

8.46 

A-74-0022-1 

15 

6.03 

5.81 

5.88 

A-7 4-0009-1 

5 

12.84 

25.53 

17.75 

A-74-0010-1 

5 

12.47 

14.02 

22.36 

A-.74-0011-1 

5 

13.38 

14.24 

23.11 

A-74-0012-1 

15 

14.34 

9.15 

7.68 

A-74-0014-1 

15 

12.79 

6.77 

3.34 

A-74-0020-1 

15 

12.51 

7.22 

6.32 

A-74-0013-1 

15 

13.88 

10.59 

12.88 

A-75-0006-1 

10 

25.15 

12.20 

21.70 

A-75-0005-1 

10 

25.73 

12.54 

40.10 

A-75-0004-1 

10 

23.66 

12.19 

30.29 

A-75-0007-1 

10 

25.21 

20.81 

84.82 

A-75-0010-1 

15 

25.61 

16.84 

14.54 

A- 75-0009-1 

15 

26.08 

8.74 

21.82 

A-75-0008-1 

15 

23.20 

9.43 

24.73 

A-74-0021-1 

15 

24.51 

13.41 

31.97 

A- 75-0019-1 

20 

25.57 

8.17 

13.07 

A-75-0018-1 

20 

26.77 

6.87 

9.63 

A-75-0017-1 

20 

29.16 

8.40 

23.39 

A-74-0018-1 

15 

8.79 

5.9i 

4.97 

A-74-0015-1 

15 

5.74 

5.12 

2.80 

A-74-0019-1 

15 

9.87 

5.76 

1.58 

A-74-0017-1 

15 

9.37 

6.86 

4.68 

A-74-0016-1 

15 

9.49 

8.45 

6.66 



zc 

3C  TTTTT 

4C  T 

5C  T 

6C  r 

7C  T 

8C 


I'  U RRRR 

W R R 
U U RRRR 

U U R R 

UtJU  R RR 


N N 
MN  N 
N N N 
N NN 
N N 


9C»* 

I nc 

II 
20 
25C 
30 
35C 
40 
50 
60 
70 
80 
85 
90 
91C 
92C 
93C 
94C 
95C 
96C 
97C 
180 
110 
IZO 
130 
140 
150 
160 
170 
180 
190 
lOllC 
1012C 
1013C 
lfll4(f 
1015c 
T016C 
1017C 
1019 
lOZOC 
1022C 
10Z4C 
1025C 
1026C 
1027C 


,o»L.»EA»  ™ri»« 

external  ruNcr 

^^X(i velocity  fpsi) 

HAX1T=50 

^ '•  iSfNG***i  fS(NG#**HAXI > Mi y i t 

mNTt«S0LUTI0N«,Xn),X(2) 

C4LL  T0RN(K,X*FK,1) 

STOPIEND 


SUBROItTtNE  FUNCT 

SURROUTINE  FlINCT  <X,FK,K> 
niMENSlON  X(2) 

00  TO  (1,2),K 

1 CONTINUE 

CALL  TURN  <K,X,FK,0) 

RETURN 

2 CONTINUE 

CALL  turn  <K,X.FK*0) 

RETURN 

END 


S U 8 R 0 U TINE  B R 0 H N 


fcSPlSMALL  number  TO  TEST  FOR  CONVERGENCE 


170 


in26C 

1040 

1050 

1060 

1070 

lOBfl 

1090 

1095 

1100 

1110 

1120 

1130 

1135 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 


xs  vector  or  initial  guesses  to  the  solution 
exteiinal  fuhct 

INTEGER  S I NO, CONV. tally, POINT (20, 20) 
niHFNSION  1SU0( 20), TEMP (20), PART (20), COE (20, 20) 
niHFNSION  X(20) 

CONV=l 
SINOsl 
NN  = N 

no  12  M--1,MAX1T 
no  1 j=i,N 

1 POlNTd  , J)  = .I 
DO  9 K=1,N 
KK=K 

IF(KK  ,GT.  1)  CALL  BCKSBSf KK , NN, X, I SUB, COE» POl NT ) 
CAl L FUNCT(X,F,KK) 

FACTORS. 001 

2 TALLY=0 

no  3 I=K.N 
ITFMP=P01NT(X, 1 ) 

HOLOsXl I TFMp) 

H=FACTOR»HOLD 

IF(ABS(H)  ,LT,  1E~7)  H=.001 
X( 1TEHP)=HOLO+H 

IF(KK  .GT.  DCALL  BCKSBS(KK,NN, X, ISUB, COE, POINT) 
CALL  FUHCT(X,FPLUS,KK) 

PART( ITeHP)s{FPLUS-F)/H 
X(ITEHP)=HOLD 

IF( ABS ( PART (lTEMP))-lE-7)26, 26,24 
24  1F(ABS(F/PART( 1TEHP))-IE20)3,3,26 
26  TALLYsTALI  Yd 

3 continue 

if(tally  .le.  n-k)  go  to  4 
factor=factor»io.o 
IF(FACT0R  .GT,  0.5)  GO  TO  14 
GO  TO  2 

4 1F(K  ,LT,  K)  go  to  5 
1F(ABS(PART(!TEMP))  ,LT.  IE-7)  GO  TO  14 
COF(K,N+l )=0 

KMAXalTEHP 
or.  TO  9 

5 KHAX=POIHT(K,K) 

DERHAXs ABS ( PARK KMAX ) ) 

KPIUS=K+1 

DO  7 I=KPLUS,N 
JSUB  = POINT(K,  I ) 

TFST=A8S(PART(3SUH)) 

IF(TEST  .IT.  DFRMAX)  GO  TO  6 

bermax=test 

POINT(KPLUS, I )=KMAX 

KHAXsJSUB 

GO  TO  7 

6 PGINT(KPLUS, I )=JSUB 

7 CONTINUE 
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lE-7)  GO  TO  14 


1540  ir<ABS(PART(KHAX)  ) .LT. 

1550  ISUB(K)=KMAX 

1560  COE(K.N+l)=fl 

1570  DO  8 JsKPLUS.N 

1580  JSUB=POINT(KPLUS, J) 


1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1870 

1880 

1890 

1891C 

1892C 

1893C 

1894C 

1895C 

18‘96C 

2000 

2020C 

2030r. 

2040C 

2050C 

2060C 

2070C 

2060C 

2090C 

2100C 

2110C 


COE(K. JSUR)=-PART( JSUB)/PART(KMAX) 

8 C0E(K,N+1)=C0E(K,N41)+PART(.ISUB)*X{ JSUB) 

9 COE(K,N  + H = (COE(K.N*l  )-n/PART(KMAX)  + X(KMAX) 
X(KMAX)=cnE(M,N+l) 

IF<N  ,GT.  DCALL  BCKSBS(N,N,X, ISUB, COE. POINT) 
IF(M, 60.1)00  TO  It 
DO  91  1=1, N 

91  IF(ABS((TFMP(1)-X(I))/X(1))  .GT.  EPS)00  TO  10 
CONV=CONV*l 
lF(C0NV-3)lt,13,l3 

10  CONV=l 

11  DO  12  1=1. N 

12  TEHP(I)=X(I) 

13  HAXIT=M 
RETURN 

14  SlNGxO 
RETURN 
END 

SUBROUTINE  RCKSBSIK. N. X, I SUB , COE . PO I NT  ) 
INTEGER  POINT(20,20) 

DIMENSION  I SUB (20 ).X( 20  ). COE (20. 20  ) 

DO  1 KHM=2,K 

KM=K-KMM*2 

K.UX»!SU8(KM-1) 

X(KMAX)=fl 
DO  2 JsKM.N 
JSU0=P01NT(KM, J) 

2 X(KMAX)=X(KMAX)*C0E(KM-1. JSUB)*X( JSUB) 

1 X(KMAX)sX(KMAX)  + COE(KH-1  ,N  + 1 ) 

RETURN 

END 


S U B R 0 II  T 


SUBROUTINE  TURN 
INPUT  DATA 


INF  TURN 
(K.X.FK.L) 


VEHICLE 

MTG 

ADIM 

RDIH 

CDIM 


GROSS  VEHICLE  HEIGHT,  NEWTONS 
niSTANCE  BETWEEN  FRONT  AXIAL  AND 
niSTANCE  between  REAL  AXIAL  AND 
HORIZONTAL  DISTANCE  BETWEEN  INSI 


VEHICLE  CG. METERS 
VEHICI E CG. METERS 
DE  METER 
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2120C 

WHEELS  AND  VEHICLE  CO 

2130C 

DniM 

HORIZONTAL  D 1 ST ANCE. BETWEEN  OUTSIDE 

METER 

21A0C 

WHEELS  AND  VEHICLE  CO 

2150C 

HDIM 

VERTICAL  DISTANCE  FROM  GROUND 

METER 

2155C 

SURFACE  TO  VEHICLE  CG 

2160C 

LDIH 

DISTANCE  POINT  OF  APPLICATION  OF  DRAWBAR  PULL  TO  CG, 

2170C 

GAMMA 

angle  BETWEEN  DRAWBAR  PULL  AND  LONG. 

AXIb  OF  vehicle 

2190C 

CTR 

RATIO  OF  CENTRIFUGAL  ACCELERATION 

- 

2200C 

(V»»2/R)  TO  GRAVITY  ACCELER A T 1 ON ( G ) 

2210C 

SPKF 

SPRING  CONSTANT  FOR  FRONT  AXLE 

NFWTONS/METER 

2220C 

SPKR 

SPRING  CONSTANT  FOR  REAR  AXLE 

nfhtons/meter 

2230C 

2240C 

TIRE 

2250C 

2260C 

SLIP 

POWERED  WHEEL  SLIP, DECIMAL 

2270C 

TIRES 

CROSS  SECTION  WIDTH  OF  UNLOADED- I NFLATFD  TIRE, METER 

2280C 

TIRED 

DIAMETER  OF  UNLOADED- INF LATED  TIRE.METFR 

2290C 

OEEL 

TIRE  deflection. DECIMAL 

2300C 

ITIR 

TIRE  CODE  - FROM  MSD  TIRE  BOOK 

- 

2310C 

TRIP 

TIRE  INFLATION  PRESSURE 

KPA 

2330C 

2340C 

2350C 

SOIL 

2360C 

2370C 

KSOIL 

n=CLAY  1=SAND 

23800 

CPR 

CONE  PENETRATION  RESISTANCE, CLAY-KPAJSAND-MPA/M 

2390C 

2400C 

PARAMETRIC 

VARIATION  FOR  SPECIFIC  VEHICLE 

2410C 

KDRI VE 

♦1=REAR  WHEEL  DRIVE-  -IsFOUR  WHEEL 

drive 

.2420C 

VEI. 

VEHICLE  SPEED,  METERS/SECOND 

2430C 

DELTA 

STEERING  ANGLE  OF  STEERED  WHEELS, DEGRFES 

2440C 

lEND 

DsHORE  data  1=LAST  DATA 

2460C 

2460C 

2470 

DIMENSION  SF(2,2),TF(2,2) 

2480 

DIMENSION  X<2) 

2481 

DIMENSION  PULL<2.2) 

2490 

DIMENSION  MT(2,2),DEFL(2,2),XNUH(2,2), ITIRGO 

(5> 

2500 

DIMENSION  TC(2,?),SC(2,2) 

2510 

REAL  LDIH 

2515 

DATA  ITIRG0/2n,29,48,62,66/ 

2520 

DATA  WTO 

, ADIM,I)DIH,CDIM,DD1M,HDIM,LDIM,GAMMA 

,SRKF, 

2530 

8 SPKR, SLIP. ITIR, TRIP. KSOIL, CPR 

254  0 

8 

, XITER 

2550 

8 /1 424  0., 1.143, 1.01 6,. 09,. 89, ,629,1.758 

,n., 25228. ,25228 

2551 

8.2,40, 

2560  A 103.5,1, 

2570  H 2.5,13292,/ 

2500  DATA  KDRIVf, VEL, DELTA, lEMD 

2590  8 /-I, 2. ,7.,!/ 

2600  5 COKTINUE 

2620  GAMHR=(3,14159*GAMHA)/10O. 

2630  DELTR=(nELTA*3.18159)/186. 
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■?640  XMASS=HTG/9,81 

2050  BETAsXd) 

2660  PSI»X(2) 

2665C  COMPUTE  PRONT  WHEEL  SIDE  SLIP  ANCLE 
2670  ALPr=DELTR-BETA-<ADIM#PSI )/VEL 

2675C  COMPUTE  REAR  WHEEL  SIDE  SLIP  ANGLE 
2680  ALPR«( (BDIM*PSI )/VEL)-BETA 

2690  CTR»VEL*PS'l/9.8l 

2700C 

2710C  PROGRAM  SEGMENT  "ROLL" 

2720C  PROGRAM  SEGMENT  TO  COMPUTE  THE  ROLL  ANGLE*  SUBSEQUENT  WEIGHT  TRANSFER* 
2730C  AND  RESULTING  NUMERIC  FOR  A 4-WHEEL  VEHICLE  NEGOTIATING  A FLAT 
2740C  HORIZONTAL  TURN  UNDER  IDEAL  STEADY  STATE  MOTION 
2750C 

2760  R0LL=(WT0/(SPKF*SPKR))*(HDIM/(CDIH**2*DDIM»*2))*CTR 

2770  ROLDEO=ROLL*180. 0/3. 14159 

2780  DEM3<ADIH«BniM)«(CDIH«DDIH) 

2790  HT(1*1)s(RDIM*DDIM*HTO/DEM)-SPKF*CDIM»ROLL 

2800  HT(1*2>*(BD1H»CDIM*MTG/DEM)*SPKF»ODIM»ROLL 

2810  HT(2.1)«(ADIH^DDIM«HTG/DEh)-SPKR*CDIM«R0LL 

2820  WT(2*2)«( ADIH*cblM*HTG/DEH)»SPKR»DDIH«ROLL 

2830  DO  82  1*1*5 

2840  IF(ITIRGO( n-ITIR)82*81*82 

2850  82  CONTINUE 

2860  STOP  "NO  MATCH  ON  ITIR" 

2870  81  GO  TO  ( 20  * 29 , 48. 62* 66  ) * I 

2880C  TIRE  6 . 0 0-1 6 * 2PR J TREAD  BUFFED  SMOOTH 

2890  20  A=15.272 

2900  B*1.04 

2910  TIREB«0.168 

2920  TIRED*0.719 

2930  GO  TO  400 

2940C  TIRE  6.00-9*4PRlTREAD  BUFFED  SMOOTH 

2950  29  A»31.8 

2960  B»l,16 

2970  TIREB*0.159 

2980  TIRED*0.516 

2990  GO  TO  400 

3000C  TIRE  7.00-16*6PR)NDCC  TREAD 

3010  48  A«52,2 

3020  6*1.0 

3030  TIREB»0.189 

3040  TIRED*0.750 

3050  GO  TO  400 

3060C  TIRE  11.00-20*12PR;TREAD  BUFFED  SMOOTH 

3070  62  A«49.5 

3080  8*0.355 

3090  TIREBsO.2854 

3100  TIRED*1.05 

3110  GO  TO  400 

3120C  TIRE  9.00-16,8PRjTREA0  BUFFED  SMOOTH 
3130  66  A*64.2 

3140  B*0.565 
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// 


3150  TIREH=n.25n 

3160  linen=0.865 

3170  <00  CONTlNue 

3180  on  500  1=1,2 

3190  DO  500  J=1 ,2 

3195C  COMPUTE  TIRE  OEFIECTION 

3200  OEfU I,J)=(n»MT(l,J)/(TRIP*A))/100. 

3210  IF(KS0IL-1  H10,42fl, 

3220  410  DEF=(1 ,-DFFL( I , J)  )*»1.5 

3230  TIRF=(J .♦(TIRER/ TIRED)  )»*0. 75 

3235C  CLAY  MORILITY  NUMBER  FOR  EACH  WHEEL 

324  0 XNUM( 1,J)  = (l00fl.»CPR»TIREB*TlRED)/(HT( I . J > *DEF* T 1 RE ) 

3250  GO  TO  500 

3255C  SAND  MOBILITY  NUMBER  FOR  EACH  WHEEL 

3260  420  XNllMf  I,.J)r(CPR»(in.»«6)»(TIREB»TIREn)»*1.5»DEFI  ( I , J )) /WT  ( I , J ) 

3270  500  CONTINUE 

3280  IF{KSOIl-n,50,50 

3290C 

3300C  50U  tire  forces  FOR  VEHICLE  IN  CLAY 

3310C 

3340  TRSL=(1.+<TIREB/TIREU) )»«.25 

3350  IFIKDRI VE-1  )21.10,10 

3360  10  CONTINUE 

3365C  FRONT  TIRE  FORCES  FOR  REAR  DRIVEN  VEHICLE  IN  CLAY 
3390  DO  15  J=l,2 

3392  XNEOW=XNUM(1 , J ) • ( 1 -2 ,26« ABS( ALPF  ) * *1 . 5 ) 

3395C  TRACTIVE  FORCE 

3400  TF(1 , J)=<12./(XNEOW**2)*.007)»WT(1, J) 

3415C  SIDE  FORCE 
3420  OTAR*ABS(ALPF) 

3422  IF<OTAR.LT,0,0073)GO  TO  12 

3424C  SIDE  FORCE  IF  ALPF  > 5 DEO 
3426  BC0N=4./nTAR»».5 

3428  SF(1 ,.!)=( 15. 4-15.4*RCON/( XNE0W-7.+BC0N) >»TF(1,J) 

3430  GO  TO  15 

3432C  SIDE  FORCE  IF  ALPF  < 5 DEG 

3434  12  RATIO  sDTAR/0.0873 

34  35  SF(1,  J > = (15.4-208. 5/<XNEOM*6.5) )*RAT10»TF(1,J) 

3440  15  CONTINUE 

3450  16  CONTINUE 

3475C  REAR  TIRE  FORCES  IN  CLAY 

3480  DO  17  J=l,2 

3490  XNE0H=XNUH(2, J)*(1.-2.26»ABS( ALPR)»»1.5) 

3500  IF(XNEOW.LE.4, )XNE0W=4. 

3510  CALL  P0HSF(TRSL,XN£0W,XNUH<2, J),SLIP,ALPR,PR,Sr(2,J),PULLC> 

3515C  TRACTIVE  FORCE 

3520  TF(2, J)=PR»WT(2, J) 

3525C  SIDE  FORCE 

3530  SF(2, J)=SF(2,J)*WT(2, J) 

3535  PULL(2, J)=PULLC*WT<2, J) 

3540  17  CONTINUE 

3550  GO  TO  300 

3560  21  CONTINUE 


3975  PULUn  , J)  = PIILLC*WT(1,J) 

3980  240  CONTINUE 

3990  on  TO  2?ll 

4000  300  CONTINUE 

4005C  COMPUTE  PRAWPAH  PULL 

4010  ZnBX»PULL(2.1)*PULL<2.2)-KDR!VE*<PULL(l.l>+PULI (li2)> 

4020  ZPBsZDBX, COS<GAMMR) 

4030  ZDBYsZnR*SIN(GAHHR) 

4035C  SUM  OF  HORIZONTAL  FORCES 

4040  0FTn0T=SF(l.l)+SF<1.2)»SF(2.1 )+SF(2.2) 

4050  S -KIIRI  VF*DELTR*{Tr  (1 ,1)4TF(1,2)  )-PSI*XMAR«:*VEI.  + Zn(-Y 

4055C  SUM  OF  moments  ABOUT  CO 

4 060  PSinOT=SF(l,t)»(ADIM4CDlM*I)ELTR)  + Sm  .2)»(  AniM-nniM*DEl  TR) 

4070  « -BOIM«{Sr (2»l)+SF<2.2n*KDRIVE*(C0IM-AniM«nELTR)#TF(l.l) 

4080  & -K0PIVE*TF(l*2)*(I)DIM4AniM»DELTR)-CDIM»TF(2.1 )4DnlM»TF(2 

4081  8.2) 

4090  8 -ZDRY*LOIM 

4100  IF(L)9fl,90, 

4110  777  PA0=VEI/PSI 

4120  HRITf(6.678)  OFLT A . VEl , PS  I . ROLOEG, R A 0 , CTR , RE T A , Al  Pf , 

4125  8 ALPR 

4130  678  FORMATdX, “STEERING  ANGLE.  OEGREES  =«.F5,1./. 

4140  8 "VEHICLE  SPEED.  METRE/SEC  =",F5.1./. 

4150  8 "YAW  VELOCITY.  RAD/SEC  G12.4./. 

4155  8 "VEHICLE  ROLL.  DEGREES  r",F5.1./. 

4160  8 "CURVATURE  RADIUS. METRE  =“.F6.1./, 

4170  8 "CENTRIFUGAL  ACC./O.  »".G12,4./, 

4180  8 "VEHICLE  SIDE  SLIP  ANGLE.  RADIAN  s"»  612. 4./. 

4190  8 "FRONT  WHEEL  SIDE  SLIP  ANGLE.  RADIAN  612.4./. 

4200  8 "REAR  WHEEL  SIDE  SLIP  ANGLE.  RADIAN  612.4.///) 

4210  WRITE(6.679) 

4220  679  F0RMAT(8X, "WHEEL". 7X. "MOBILITY". 4X. "TRACTIVE". 4X. 

4230  8 "SIDE". /.ex. "LOAD. N",6X. "NUMBER". 6X, "FORCE. N".5X. 

4240  8 "FORCE. N") 

4250  PRINT  1000.((II.JJ.HT(II.JJ),XNUM(II.JJ).TF(II.JJ). 

4260  8 Snn..lJ).JJ«1.2).II*1.2) 

4270  1000  F0RMAT(2H  < . Il.lH.  . I1.1H)’.4612.5) 


4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350C 

4360C 

4370C 

4380C 

4390C 

4400C 

4410C 

4415C 

4416C 


IF(K.NE.l.aR.K.NE.2)RETURN 
90  GO  TO  (1.2).K 

1 FKxBETOOT 
RETURN 

2 FKxPSIDOT 
RETURN 
END 


POWSF  SUBROUTINE  TO  COMPUTE  SIDE  FORCE  COEFFICIENT 
(S/W)  FOR  POWERED  WHEELS  OPERATING  IN  CLAY 


" -ar--- 


3575C  TRONT  TIRE  FORCES  FOR  ALL  WHEEL  DRIVE  VEHICLE  IN  CLAY 
3590  DO  22  J=l,2 

360  0 XNFOH=XNUH(l , J)*(1.-2.26*A0S( ALPF»»»1 15) 

3610  IFIXNEOM.I E.4.  )XNEOHs4. 

3620  CALI  POMSFITRSL. XNEOH.XNUHd, J),SLIP. ALPF,PF,SF(1, J),PULLC) 

3625C  TRACTIVE  FORCE 

3630  TF<1  . J)=1>F.HT(1 , J) 

3635C  SIDE  FORCF 

3640  SF(ti.n  = sm.  n»HT(i,.j) 

3645  Pin  I (1  , J)=Rilll  C«WT(1  . J) 

3650  22  CONTINUE 

3660  GO  TO  36 

3670  50  CONTINUE 

3660C 

3690C  SOIL  TIRE  FORCES  FOR  VEHICLE  IN  SAND 
3700C 

3740  IFIKDRI  VE-1  )?5n.2n0, 

3750  200  CONTINUE 

3755C  FRONT  DRIVE  FORCES  FOR  REAR  DRIVEN  VEHICLE  IN  SARD 
3760  DO  210  J=1 ,2 

3765C  TRACTIVE  FORCE 

3770  TF(l, J)=( .0t5*.837(XNUh(l. J)-2. > )»HT(l, J) 

3775  DTARxAHSIALPF) 

3777  IFLDTAR.LT, 0.0073)00  TO  195 

3779C  SIDE  FORCE  IF  ALPF  > 5 DEG 

3780  SF<1. J)  = (1 .275»DTAR»»1 .23+.83-46./(XNUM(l, JI+55.4 ) )*HT(1, J) 

3782  00  TO  190 

3783C  SIDE  FORCE  IF  ALPF  < 5 DEG 

3784  195  RATIO=DTAR/n.0873 

3706  SF(1, J)=RATIO*(.09-46./(XNUM(l, J)*55.4) )*HT(1».I) 

3780  198  CONTINUE 

3790  210  CONTINUE 

3800  220  CONTINUE 

3815C  REAR  TIRE  FORCES  IN  SAND 

3830  DO  230  J»1 ,2 

3840  DHMALPsABS(ALPR) 

3850  CALL  PHSAN ( XNUM( 2, J ) , SL I P. DIIMALP. TC ( 2» J ) , SC ( 2* J ) . PULLC ) 

3855C  TRACTIVE  FORCE 

3860  TF(2, J)=T0(2, J)*HT(2, J) 

3865C  SIDE  FORCE 

3870  SF(2, J)=SC(2. J)*HT{2,J) 

3875  PULI  (2.  J)  = PULI.C»HT(2*  J) 

3880  230  CONTINUE 

3890  GO  TO  300 

3900  250  CONTINUE 

3915C  FRONT 'TIRE  FORCES  FOR  ALL  WHEEL  DRIVE  IN  SAND 
3930  DO  240  J = 1 ,2 

3940  DUMARxABSCALPF) 

3950  CALL  PWSANI XNUMd , J ) , SL I P. DUMAR, TC ( 1 , J > . SC ( 1 , J ) , PULLC ) 

3955C  TRACTIVE  FORCE 

3960  TF(1,J)=TC(1, J)»WT<1,J) 

3965C  SIDE  FORCE 

3970  SF(t, J)=SC(1 , J)*WT<1, J) 
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4420  ••>.‘SWB?>OUHNE  POWSF (TRSL.XNEOK^XMUH, SLIP, ALP* TPOW.SFP.PULLO 

4430  ‘ -SSPaiai./IXNEON^^a.SMt.OOS 

4440  SLr»SLIP/SSP 

4450  PULLC«.5*AL0010(SLD 

4452  0SPa(12./XNE0N«*2. )+.fl07 

4453  0=TRSL«PUl.Ln  + OSP 

4454  IF(PULLC.(?T.Q)00  TO  10 

4455  TPOW»0 

4456  GO  TO  20 

4457  10  TPOH«PULLC 

445fl  20  continue 

4460  SIX  <XNUH44. )«o2/16.«0.00B1>**0.5*(3.37*ARS(ALP) 

4470  A -4.24*ABS(ALP)*»2)»(l./.46) 

4480  SL0PE*1.7189*ABS(ALP) 

4490  SFP=SI-TPOW*SLOPE 

4500  RETURN 

4510  end 

4520C  

4530C 

4540C  SUBROUTINE  PMSAN 
45S0C 

4560C  •••*•••••• 

4570C 

4575C  PNSAN  SUBROUTINE  TO  COHPUTE  SIDE  FORCE  COEFPICIENT(S/M) 
4576C  and  TRACTTVF  FORCE  COEFFICIENT  FOR  POWERED  WHEELS 

4577C  OPERATION  IN  SAND 

4580  SUBROUTINE  PWSANIXNUN, SLIP, ALP, TCS,SCS.PULL'C) 

4590  A«.69-.01/SLIP-1.42*ALP**(.6*3.1*SLIP) 

4600  B>10.8*16.*ALP«*.8 

4610  C«2.23/SLIP**<l./3. )^<15./SLIP)*ALP**3. 

4620  PULLC>(A-A«e/(XNUM>C«B)) 

4621  BH«4,71+1.72/SLIP 

4622  0«.66-.66*BH/(XNUN«10.«BN) 

4623  IF(PULLC.6T*0)00  TO  10 

4624  TCS«0 

4625  60  TO  20 

4626  10  TCS«PULLC 

4627  20  CONTINUE 

4630  Ea2.34.03*XNUN 

4640  F«2.4».065*XNUH 

4650  SZPU«E*ALP-F*ALP**2. 

4660  SCSaSZPU-(1.65»ALP>«TCS 

4670  RETURN 

4680  END 


In  accordance  with  HR  70-2-3,  paragraph  6c(l)(h), 
dated  Fehruazy  1973,  a facsisile  catalog  card 
In  Library  of  Congraae  fonaat  Is  reproduced  below. 
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